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ABSTRACT 

The zCOSMOS SINFONI project is aimed at studying the physical and kinematical properties of a 
sample of massive z ~ 1.4-2.5 star-forming galaxies, through SINFONI near-IR integral field spec- 
troscopy (IFS), combined with the multi-wavelength information from the zCOSMOS (COSMOS) 
survey. The project is based on 1 hour of natural-seeing observations per target, and Adaptive Optics 
(AO) follow-up for a major part of the sample, which includes 30 galaxies selected from the zCOS- 
MOS/VIMOS spectroscopic survey. 

This first paper presents the sample selection, and the global physical characterization of the target 
galaxies from multicolor photometry, i.e., star formation rate (SFR), stellar mass, age, etc. The Ha 
integrated properties such as, flux, velocity dispersion, and size, are derived from the natural- seeing 
observations, while the follow up AO observations will be presented in the next paper of this series. 
Our sample appears to be well representative of star- forming galaxies at z ~ 2, covering a wide range 
in mass and SFR. The Ha integrated properties of the 25 Ha detected galaxies are similar to those 
of other IFS samples at the same redshifts. Good agreement is found among the SFRs derived from 
Ha luminosity and other diagnostic methods, provided the extinction affecting the Ha luminosity is 
about twice that affecting the continuum. A preliminary kinematic analysis, based on the maximum 
observed velocity difference across the source, and on the integrated velocity dispersion, indicates that 
the sample splits nearly 50-50 into rotation-dominated and velocity dispersion-dominated galaxies, in 
good agreement with previous surveys. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: kinematics and dynamics 
— infrared: galaxies 
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1. INTRODUCTION 

For some time, the high star formation rates (SFR) and 
the ensuing stellar mass growth of massive galaxies at 
early cosmic epochs have been widely attributed to vi- 
olent major mergers. Motivated by a number of recent 
observational findings, our picture is now changing to- 
wards one in which smoother yet efficient modes of gas 
accretion fueling intense star formation play an impor- 
tant, if not dominant, role at z ~ 2. In fact, spatially- 
and spectrally-resolved integral-field spectroscopy has 
revealed large rotating disks at z ~ 2 with SFRs as high 
as ~ 100 M yr~ l , or more, but without any sign of ongo- 
ing major merging ([Forster S chreiber et al. 2006; Forster 
Schreiber|[2009l |Cresci et al.||2009t |Genzel et al.| |2006 , 



2008) |2011 UShapiro et al.|2008| ). In parallel, multiwave- 
length galaxy surveys have shown that the SFR correlates 
tightly with stellar mass and is nearly proportional to 
it, with a small scatter (~ 0.3 dex in logarithmic units), 
and that the overall SFR-M^ relation steadily declines 
from z ~ 2.5 to z ~ ([Daddi et al.|2007||Elbaz et al.|2Q07 



Noeske et al. 2007 



|2Q09| [SantTni et al. 
somewhat flatter S 



IPannella et al.||2009[ |Dunne et aL] 
2009[|Damen et al.||2009|), although 



R-M^ relations have also been de- 
rived (Rodighi ero et al.||2010[ |Karim et aL|2011| |Wuyts 



et al.|[2UTTT T Moreover, recent high spatial resolution 
observations of CO emission from "normal" (i.e., non- 
major merging) massive star- forming galaxies at z ~ 1 -2 
have uncovered large molecular gas mass fractions, three 
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to ten times higher than in local massive disks, and indi- 
cate that the star formation efficiency (i.e., the surface star 



not evolve strongly with cosmic time ( 


Daddi et al.|2010| 


|Tacconi et al.|2010||Genzel et al.|2010 





Thus, the majority of star- forming galaxies appears to 
be continuously fed by gas, promoting and maintaining 
star formation, rather than occasionally bursting as a re- 
sult of a (major) merger. The regularity and simplic- 
ity of the SFR(M^,f) relation all the way to z ~ (e.g ., 
Peng et al.|2010UBouche et al |2010[|Dutton et al.|2 010), 



the mass-metallicity relation ( |Erb et al. 2006a]), and the 
tight fundamental metallicity relation (FMR), between 
Af*, metallicity, and SFR (Man nucci et al.||2&T0] ), actu- 
ally suggest that this gas accretion process is the dom- 
inant mechanism driving the mass growth of galaxies 
over a wide range of redshifts. This new empirical evi- 
dence matches remarkably well with recent hydrodynam- 
ical simulations in which massive galaxies acquire a large 
fraction of their baryonic mass via quasi steady, narrow 
cold flows or streams that penetrate effectively through 
the s hock-heated media of massive dark matter halos 
(e.g., IBinney|[2004l |Keres et al.|[2005| |Dekel and Birn 



|boim|2006[ |2008; Deke l et al.|2009| ). These cold streams 
consist of smoothly flowing material (whose clumps may 
be seen as minor merger events) that can sustain elevated 
SFRs over much longer timescales than major mergers. 
Under these conditions, the angular momentum is largely 
preserved as matter is accreted, early (thick) disks can 
survive and be replenished, and internal dynamical pro- 
cesses can drive secular evolution of disks and forma- 
tion of bulges/spheroids at high z. Indeed, fragmenta- 
tion in such turbulent, gas-rich disks tend to form mas- 
sive self-gravitating clumps that could migrate inwards 
and coalesce to form a young bulge in ~ 1 - 2 Gyr, as 
argued both on observationa l (e.g. |Genzel et al.| 2008 , 
|2011[|Elm egreen et al. 2009 ]) and on theoretical grounds 
(e.g. Noguchi 1999 ; Bournaud and Combes 2002 ; Immeli 
et al.|2004a|bhCarollo et al.|2007t|Bo urnaud et al. 2007; 
Elmegree n et al.|20 08 ; Deke l et al.|20 09 ; Ceverino et al.| 



2010p . In contrast, other studies suggest that giant clumps 
may not survive long enough to migrate to the galaxy cen- 
ter since they could be rapidly disrupted by vigorou s stel- 
lar feedback ( [Murray et al.|2010||Gehel et al.|2010] ). 

Spatially -resolved gas kinematics of high redshift ob- 
jects has provided key insights into the mechanisms driv- 
ing early galaxy evolution in the above framework. The 
"SINS" project, carried out with the SINFONI integral 
field spectrograph at ESO VLT, was the first and largest 
survey with full 2D mapping of the spatial distribution 
and kinematics of the Ha line emission for 62 massive 
z ~ 1.3-2.6 star-forming galaxies, on typical resolved 
scales of ~ 4-5 kpc and down to ~1.5 kpc for a sub- 



set ob served with adaptive optics (AO) (Forste r Schreiber 
|2009[ hereafter FS09, and references therein). About 1/3 
of the SINS galaxies exhibit clear kinematic signatures 
of ordered rotation in a disk, despite generally irregular 
Ha morphologies. Key properties of these early disks 
are large intrinsic velocity dispersions and high gas mass 
fractions of ~ 30% or more. Concurring evidence has 
been found fr om other dynamical and/or morphological 
studies (e.g., [W right et al.|2007[ [20091 [El megreen et al. 
20071 |Bournaud et al.|2008]|van S tarkenburg et al. 2008; 
Law et al.||2007[ |2009[ |Stark et al.||2008[ |Epinat et al.| 



2009 1 |Jones et al.|20T0l |Gnerucci et al.|2011| ). 

While the results from SINS and other recent studies 
have provided new and significant insights into galaxy 
formation/evolution, the picture remains very incom- 
plete. Indeed, another key outcome of SINS and other 
IFS studies is the diversity of dynamical properties 
among the sample galaxies. Along with large massive 
disks, these studies have revealed a population of com- 
pact gas rich systems, frequent at lower masses, that ap- 
pear dominated by random moti ons at the current resolu- 
tion of 8m class telescopes ( [Law et al.|2QQ7] [2009] FS09), 
as well as a number of merging/interacting systems. The 
nature of the dispersion dominated galaxies, and the re- 
lationship between all three dynamical classes of high-z 
galaxies are still unclear. The present new series of SIN- 
FONI observations have been conceived to address these, 
and other related open issues, such as: why these objects 
do follow the same SFR(M*,0 relation, in spite of their 
morphological and dynamical diversities? and how they 
relate to the cold-stream paradigm? In an even broader 
perspective, with these SINFONI observations we would 
like to start investigating issues such as the formation 
of galactic bulges, star formation quenching in massive 
galaxies, and ultimately how the diversity and complex- 
ity of the physical processes at work may result in the 
remarkable simplicity in the evol ution of galaxy po pula- 
tions as extensively illustrated by Pen g et aL] ( |2010| ). 

To make progress, it is now important to expand the 
currently still limited sample of galaxies at z ~ 2 with 
high S/N information on the kinematics and distribution 
of star formation resolved on ~ 1 kpc scales. High sensi- 
tivity and high spatial resolution are necessary to address 
the ultimate science goals outlined above. More specif- 
ically, these goals require to study separately structural 
components within galaxies (disk/bulge regions, massive 
star-forming clumps), resolve merging units and more 
compact lower-mass systems, and detect kpc-scale per- 
turbations in the velocity fields and dispersion maps as 
well as faint emission line profile components signal- 
ing dynamical interactions or processes such as radial in- 
flows/outflows. 

This can be achieved with AO-assisted near-IR inte- 
gral field spectroscopy with spectral resolution of R = 
A A/ A ~ 3000-5000 but it remains prohibitive in terms 
of telescope time to collect such data for samples of 
galaxies of size comparable to that of current spectro- 
scopic surveys. 

Our aim is therefore to observe with SINFONI and AO 
a controlled set of galaxies at z ~ 2, i.e., around the peak 
of the cosmic SFR density (e.g. |Hopkins and Beacom 



2006 ; |Le Borgne et al.|2 009 ) and AGN activity ( [Fan et al. 
|200l[ |Croom et al.||2004| ), with well-characterized pho- 



tometric, spectral, and morphological properties. This 
will allow a detailed study of their kinematics in con- 
junction with their stellar mass, star formation, dust ex- 
tinction, and structural properties. This set of benchmark 
z ~ 2 star-forming galaxies will help to understand the 
processes at play at the heyday of massive galaxy for- 
mation in the much broader sample of galaxies in cur- 
rent spectroscopic and photometric surveys. In this per- 
spective, we have collected a set of 30 1.4 < z < 2.5 
massive s tar-forming galaxi es from the zCOSMOS-Deep 
database (Lilly et al. 2007, hereafter L07) spanning the 
mass range 3 x 10 y < (Af*/Af ) <2x 10 11 and prob- 
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ing the "main-sequence" of star-forming galaxies over 
this mass range. zCOSMOS-Deep is the most appro- 
priate survey in terms of the combination of field size 
and spectroscopic completeness at z ~ 2. It thus pro- 
vides the best opportunity to cull our benchmark sample, 
which needs to satisfy the stringent criteria for feasible 
AO-assisted SINFONI observations: accurate redshifts at 
z ~ 1.4-2.5, and the proximity of a bright star for AO 
correction. In addition, the extensive multi wavelength 
observations available for the COSMOS field from the 
X-ray to the radio regimes provide rich complementary 
data sets on each target, as well as for other ~ 10,000 star 
forming galaxies in the same redshift interval. 

Our SINFONI programs includes two parts: (i) the 
natural- seeing " pre -imaging" observations, carried out 
without adaptive optics correction (no-AO mode) for a 
total sample of 30 zCOSMOS star-forming galaxies at 
1.4 < z < 2.5, typically with one hour of total integration 
time per object, and (ii) the AO follow-up observations 
for a subset of about 20 sources, selected based on the 
strength of the Ha line emission in the pre-imaging data, 
and with integration times in the range ~ 4-10 hours per 
target. Combined with AO data of another eight galax- 
ies drawn from the SINS survey, the total benchmark AO 
sample will reach nearly 30 objects for an increase by fac- 
tors of at least two over currently published AO samples 
at z ~ 2. In this paper we present the general properties 
of the full sample, along with the first results of our SIN- 
FONI+VLT pre-imaging observations consisting of one 
hour integrations under natural seeing conditions. The 
follow-up AO-assisted SINFONI observations are now 
well under way, and will be presented in a series of forth- 
coming papers. The present paper is primarily meant to 
fully characterize the target galaxies, and compare them 
to the general galaxy population at the same redshift. 
This new sample of 30 objects with natural- seeing SIN- 
FONI data expands significantly on similar existing sam- 
ples at z ~ 2 and constitutes the second largest one after 
the SINS Ha survey of 62 objects (FS09). The results 
from the seeing-limited SINFONI data discussed in this 
paper are therefore valuable in their own right, notably in 
view of future seeing-limited surveys of larger samples 
with the near-IR multi-integral field spectrograph KMOS 
at the VLT. 

The paper is organized as follows. Section [^de- 
scribes the sample selection criteria and Sectionp] il- 
lustrates the properties of each galaxy in the sample as 
derived by their photometric spectral energy distribu- 
tion (SED), and compare them to the properties of other 
similar galaxy samples from the literature. Section [4] 
focuses on the SINFONI observations, data reduction 
and the Ha line integrated properties. In Section [5] we 
compare the SFR derived from the Ha luminosity with 
those obtained through other SFR indicators. Finally, 
in Section [6] we perform a classification into rotation-, 
and dispersion-dominated sources, based on the seeing- 
limited pre-imaging observations. Section [7] is the sum- 
mary and conclusions. Throughout this paper we adopt 
a cosmology with Ho=71 kms -1 Mpc -1 , £7m=0.27, and 
fi A =0.73. 



2. TARGET SELECTION: THE ZC-SINF SAMPLE 

The zCOSMOS-SINFONI sample (hereafter zC-SINF) 
was extracted from the spectroscopic zCOSMOS-Deep 



surve\ 



? y coll ected with VIMOS at the VLT (L07; |DHy| 
et al. 2009; and in prep.), covering t he central 1 deg 2 of 



the whole 2 deg 2 COSMOS Survey ( |Scoyille et al.|2007 ). 
With about 10,000 spectra having been obtained, and reli- 
able redshifts derived for over 50% of them, zCOSMOS- 
Deep is the largest spectroscopic redshift survey of 1 .4 < 
Z < 3 galaxies. It is part of the larger zCOSMOS pro- 
gram, including also the 'zCOSMOS-Bright' component, 
i.e. ~20,000 galaxy spectra covering the whole COS- 
MOS field area (1.7 deg 2 effective area), a nd the redshift 
range 0.1 < z < 1.2 (L07 |Lilly et al.|2 009, and in prep.). 

For the early zCOSMOS-Deep observations, under- 
taken in 2005, targets were selected using a relatively 
shallow ^T -band imaging and catalog (K AB < 21.85, |Ca^| 
|pak et al.||2007) . However, the bulk of the zCOSMOS- 
Deep targets were observed in 2007 (see L07), by tak- 
ing advantage of the deeper COSMOS Wircam+CFHT 
AT-ba nd imaging, and of the rel ative ^T-band selected cat- 
alog (Mc Cracken et al.||2Q10] ). This catalog is > 70% 
complete down to K AB = 23 lor extended sources (> 50% 
down to K AB = 23.5). The pre-selection of zCOSMOS- 
Deep spectroscopic targets was based on two-color crite- 
ria, both v ery efficient in sel ecting 1.4 < z < 2.5 galaxies: 
the BzK paddi et al.pOO?]), and th e UGR BM/BX opti- 
cal color criteria ( Steidel et al. 2004). The combination of 
both criteria was initially necessary to ensure a sufficient 
surface density of suitable targets, to take full advantage 
of the VIMOS instrument multiplex. The BzK criterion, 
coupled with a deep AT-band selection (i.e., K < 23.5, 
for zCOSMOS-Deep), enables one to select both passive 
galaxies (pBzK) and star-forming galaxies (sBzK, across 
a range of dust obscuration) at 1.4 < z < 2.5. The most 
massive (^-brightest) galaxies are largely missed b y op- 
tical (UV rest-frame) color selection criteria (seejDaddi 
|et al. 12004) for more details). The sBzK-selected objects, 
the relevant star- forming population in the context of this 
work, have on average higher mass, reddening, and SFR 
compared to the U GR- selected galaxies at the same red - 
shift (z ~ 2) (cf. |Daddi et al.||2004[ |Reddy et al.| |2005 ). 
Whereas the BzK criterion with deep K-bmid data picks 
also the vast majority of the BM/BX objects, it could 
miss some younger and less massive galaxies in the red- 
shift range of interest, and therefore the UGR BX and 
BM selection criteria were also used to recover this mi- 
nor part of the population. In zCOSMOS -Deep 5-band 
magnitude cuts were also applied to ensure an adequate 
S/N ratio in the optical continuum of VIMOS spectra, i.e., 
B AB < 24.75, and B AB < 25.25 for sources with photomet- 
ric redshift z p hot < 2, and z p hot > 2, respectively. This dif- 
ferent depth was dictated by the need of measuring red- 
shifts from relatively low S/N spectra which is harder in 
the lower redshift interval when Lyman-a is not included 
in the VIMOS spectral range. Moreover, the zC-SINF tar- 
gets were selected to match the following requirements: 

1. Proximity of a star suitable for natural guide star 
(NGS) adaptive optics correction (with g AB < 17 
mag, and within 30" of the galaxy). A total of 
622 of the photometric candidate objects satisfied 
this criterion, forming the basis for culling the SIN- 
FONI targets according to the next three selection 
steps. 

2. Spectroscopic redshift reliability: only spectra with 
the best L07 Confidence Classes, e.g., 4, 3, and 
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2 were used. The three classes indicate, respec- 
tively, "completely secure redshifts" (i.e., based on 
unambiguous multiple spectral features of the ex- 
pected relative strength that leave no room for any 
doubt about the redshift), "very secure redshifts" 
(for which the classifiers recognize a remote pos- 
sibility for error, e.g., because supporting features 
are in a noisy part of the spectrum), and "less se- 
cure redshifts" (i.e., for which the claimed redshift 
is the most likely, but a not negligible possibility 
remains that the redshift is not correct). A subset 
of 249 galaxies out of the 622 candidates close to a 
suitable NGS reference star met this requirement. 

3. Ha line "observability:" we then considered only 
the sources with redshifted Ha line falling either 
in the SINFONI H- or 7^-band, within spectral 
regions of high atmospheric transmission, and at 
least 400 kms~ l away from OH airglow lines with 
typical intensity > 10% of that of the brightest one 
in either H or K band as appropriate. 

This criterion constrains the target redshifts either 
in the range z ~ 1.3 - 1.7, or z ~ 2-2.5, respec- 
tively, as shown in Figure [Tj where the spectro- 
scopic redshift distribution of the zC-SINF sample 
(red histogram) is compared with the photometric 
redshift distribution of the total sample of COS- 
MOS BzK gal axies at 1.4 < z < 2.5 (photometric 
redshifts from Ilb ert et aL]|2009| ). Together with 
the OH avoidance requirement for Ha, this led to a 
sample of 65 candidates. 

4. Based on modeling of the opti cal to mid-IR spec- 
tral energy distribution (using Bruzual and Char- 

| lot||2 003 with constant SFR, see next Section) we 
selected targets so as to probe the widest possible 
range in stellar mass (~ 3 x 10 9 -2 x lO n M ) and 
SFR (~ 10-300M©/yr). The lower SFR limit was 
meant to ensure Ha detection in ~ 1 hr integra- 
tion time with SINFONI i n seeing-limited mode 
(see, e.g., Forster Schreiber 2009, hereafter FS09), 
and the resulting sample of candidates included 62 
galaxies. 1 

The final zC-SINF sample we observed with SINFONI 
is composed of 30 zCOSMOS-Deep sources. Among the 
62 candidates that satisfy all criteria listed above, they 
represent the 30 best ones in terms of combined optical 
redshift reliability, Ha line observability, and coverage 
in and SFR parameter space. Four targets were ob- 
served in earlier campaigns as part of MPE SINFONI 
guaranteed time observations (GTO) in the spring of 2007 
(P79), and used as a 'pilot sample' for the project (see Ta- 
ble [TJ. These four objects are al so part of the 'SIN S Ha 
sample' presented by FS09 and |Cresci et al. ( 2009). Al- 
though we did not apply any selection criterion based on 
the galaxy size, relatively extended, low surface bright- 
ness sources are likely to suffer a negative bias, when se- 
lected as zCOSMOS targets on the base of magnitudes, 
as for them getting redshifts from VIMOS data was likely 
harder than for more compact sources. 

1 The impact of the (a priori unknown) Ha size or surface brightness 
distribution on source detectability at a given integrated line flux or SFR 
is discussed in § |4.3.3| 



0.5 



0.4 



0.3 



0.2 




o.o h m , i 



Fig. 1 . — The redshift distribution of the sample of 30 zC-SINF 
targets (red dashed histogram, right scale), among which 25 
have Ha detection with SINFONI (red filled area), is compared 
with the photometric redshift distribution of all the ~ 30, 000 
BzK- selected galaxies in COSMOS (normalized black dashed 
line histogram, left scale) . The latter histogram was built by 
ma tching the high-z Bz KpVlcCracken et al. ( 2010 ) catalog with 
the |Ilbert et aL] p009| ) photometric redshifts! The observed 
peaks around z ~ 1.5 and z ~ 2-2.5 in the red histogram are 
produced by the constraint of having the Ha line within the 
SINFONI H- or ^T-band, respectively, combined with the other 
selection criteria, as discussed at the end of Section[2] 

As shown in Figure [T] 25/30 objects are detected in 
Ha (red histogram). The remaining objects (dashed area 
in Figure [T]) are undetected possibly due to higher noise 
levels associated with the proximity in wavelength of 
night sky lines 2 , Ha fluxes lower than expected from the 
spectral energy distributions modeling, surface bright- 
ness limitations for more extended sources, or a com- 
bination of these effects. The redshift distribution of 
zC-SINF targets, shown in Figure [T] is much different 
from that of the general COSMOS population (dashed- 
line histogram), due to a combination of biases coming 
from the selection criteria described above. First, the re- 
quirements to detect Ha in the SINFONI H- or ^T-band, 
coupled with the requirement to avoid atmospheric con- 
tamination to Ha excludes galaxies in the 1.6 < z < 1.9 
redshift range. Moreover, the peak observed at z ~ 2.1 
reflects the peak of the spectroscopic redshift distribu- 
tion of zCOSMOS-Deep galaxies (Lilly et al in prep.), 
resulting from the convolution of the actual redshift dis- 
tribution of zCOSMOS-Deep targets with the success 
rate in redshift determination. Indeed, the most visible 
consequence of this issue, is the deficit of galaxies at 
z = 1.4- 1.6 in the final zC-SINF sample (only including 
two of them). Nevertheless, considering that in the SINS 
sample there is a fair number of galaxies in the redshift 
range 1.4 < z < 1.6, these two galaxies should not be re- 

2 Despite our OH avoidance criterion, it is possible that for some 
cases, weaker OH lines affected the noise levels sufficiently to hamper 
detection of a faint Ho; emission line. It is also possible that the optical 
redshift differs from the Ha redshift by a few 100 km/s (as for some of 
the detected objects, see Table 5) so that the Ho; line falls closer than 
expected to a bright OH line, preventing detection. 
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garded "in isolation", but as complementing and expand- 
ing the sample of galaxies with resolved Ha kinematics 
at these redshifts. 

3. PROPERTIES OF THE TARGETS COMPARED TO THE 
BROADER POPULATION OF STAR-FORMING GALAXIES AT 

Z~2 

In this section we present the multi-band photometric 
data, and the general properties of the zC-SINF sample, 
as derived from the Spectral Energy Distribution (SED) 
analysis. We test the reliability of masses and SFRs de- 
rived through SED fitting against the assumptions on stel- 
lar population models, star formation histories, and pho- 
tometric catalogs. Then, we compare the zC-SINF galaxy 
properties with those of the whole sample of z ~ 2 COS- 
MOS galaxies to establish whether our sample is rep- 
resentative of the star-forming galaxy population in the 
same redshift range. 

3.1. Multi-band photometry 

We used the multi-band photometry from the deep 
WirCam+CFHT (Canada -France-Hawaii Telescop e) K 
band selected catalog of McCracke n et al.j (|2010|), in 



tegrated wit h the IR/Mid-IR photom etry of jllbert et al. 
(2009| and |Le Floc'h et~aLl ([2009] ). While the pub- 
lished [McCracken et al. ~ pQ10| ) catalog provides mag- 
nitudes for few bands (i.e., Bj, i + , z + , /, and K, see 
that paper for more details), in this work we used an 
updated version of the catalog (E. Daddi, and H.J. Mc- 
Cracken private communication), complemented with the 
WirCam+CFHT //-band, and four more optical bands, 
i.e., the Mega-Prime+CFHT w*-band, and the Suprime- 
Cam+ Subaru Vj-, and r + -band data (Capak et al. 
2007| ). The u*, g + , V/, r + , and H photom etry was de~ 
rived using the same procedure described in |McCracken| 



et aLl pOTO] ) for the other bands. Magnitudes were mea 
sured in 2 diameter aperture, and corrected to total by 
applying an offset derived for point-like sources within a 
6" diameter aperture. For the four new optical bands, the 
following corrections were applied: 



u* ot = u* -0.212 
g ; ot = g + - 0.245 

v,, M = y,-o.i6 

r+ ( = r + -0.125 

The //-band aperture correction depends sensitively on 
the position in the COSMOS field, varying from -0.145 
to -0.215 magnitudes across the field (E. Daddi private 
communication) . 

For the other optical and near-IR bands we ref er to the 
corrections reported in McCrac ken et al.| ( |2010| ). How- 
ever, the above corrections provide the total magnitude 
for point-like sources, while extended objects need a fur- 
ther correction obtained by comparing the ZT-band aper- 
ture total mag nitude with the MA G_AUTO provided by 
SExtractor ( [Berlin and Arnouts1|1996| ), i.e., the magni- 
tude measured in a Kron-like elliptical aperture (includ- 
ing more than 90% of the flux), and generally called 
'auto-offset'. The 'auto-offsets' derived in ZT-band for 
each galaxy, were equally applied to all the optical and 
near-IR bands, and, for our sample, they typically range 
from -0.05 to -0.35 magnitudes. As mentioned above, the 
catalog was cross-matched with the four Spitzer+IRAC 



channels of Ilbert et al. (2009), covering from 3.6 fim to 
8.0 //m, and with the 2 4 /mi Spitzer+MIPS deep catalog 
( |Le Floc'h eTaTl|2009] ). The final multi-band photome- 
try from optical to mid-IR for all the zC-SINF sample is 
reported in Table [2] and [3] 

3.2. SED fitting analysis 
3.2.1. Assumptions and procedures 

The optical-to-IRAC photometric data reported in Ta- 
bles]^ and [3] were used to perform SED fitting analysis to 
estimate galaxy stellar mass, SFR, dust extinction (Ay), 
and age (i.e., the time elapsed since the beginning of star 
formation). In a few cases the sources were blended in 
some of the IRAC bands (due to the poorer angular res- 
olution compared to the other bands), and such IRAC 
data were excluded from the fit. We used the HyperZ- 
mass softw are, i.e., a modified ver sion of the public Hy- 
perZ code (jBolzonella et al.J2000| ), which fits photomet- 
ric data points with synthetic stellar population models, 
and picks the best fit parameters by minimizing the x 2 
between observed and model fluxes. The stellar mass is 
obtained by integrating the star formation rate over the 
galaxy age, and correcting for the mass loss in the course 
of stellar evolution ( ~ 40% for a galaxy o f ~ 1 Gyr and 
a Chabrier IMF) (see |Pozzetti et al.||2007] ). The best fit 
results strongly depend on the assumptions concerning 
metallicity, initial mass function (IMF), stellar population 
models, extinction law and star formation history (SFH). 
In order to minimize the number of free parameters in 
this work we used onl y models built with solar metal- 
licity, and adopted the lChabrier| ( |2003| ) IMF. We fixed 
the redshift t o the spectroscopic values, used thelCalzetti 
et al. (2000]) la w to account f or dust extinction, and the 



prescription of |Madaul ( |1995| ) to treat the Lyman-a for- 
t. The dust < 



est flux decrement. The d ust extinction param eter (i.e., 
A v = E(B - V) x 4.05, cf. |Calzetti et alT1[2000] ) was al- 
lowed to range from Ay = to Ay = 3 in steps of 0.1. 



As recently discussed in Mara ston et al.| ( 20T0| here- 
after MA10), the assumption on the SFH is crucial to 
derive astrophysically plausible best fit parameters, and 
should be chosen based on the expected properties of the 
galaxy sample. In particular, MA 10 argued that models 
with constant SFR (CSFR), and the so-called "r-models" 
characterized by a SFR exponentially decreasing with 
time as e~ t ' T , do not give a realistic representation of z ~ 2 
star-forming galaxies. When age is left as a free parame- 
ter both models lead to very low ages (< 0.1-0.3 Gyr), 
which would imply a formation redshift almost identical 
to the redshift at which the galaxies are observed. As ex- 
tensively discussed in MA 10, such unrealistically young 
ages are the consequence of the outshining effect by the 
most recently formed stellar populations, and should not 
be considered real. In an attempt to circumvent this diffi- 
culty, a lower limit (of ~0. 1 Gyr) is often imposed to the 
galaxy age, with the result that many (most) galaxies are 
found to cluster at such imposed minimum age. This is 
because in such fits the current SFR and Ay are dictated 
by the rest-frame UV part of the spectrum, but within 
the assumed SFH such a SFR would largely overproduce 
mass (then violating the near-IR flux constraints) if not 
for very short ages. Such unrealistically short ages suffice 
to demonstrate that the average past SFR must have been 
much lower than the current values, hence SFR (on aver- 
age) must has been secularly increasing with time. Thus, 
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the best-fit ages derived using either CSFR or r-models 
should rather be interpreted as those of stars providing the 
bulk of the light in those bands for which the photometry 
is affected by the smallest errors (as it is built-in in the x 2 
procedure), i.e. in the optical, which corresponds to the 
rest- frame UV. 

Moreover, "r-models" automatically assume that all 
galaxies are caught at the minimum of their SFR, while 
it appears more plausible that z ~ 2 star-forming galax- 
ies are close to their SFR peak, given that this is the 
cosmic epoch of maximum star formation activity and 
given the function al form of the SFR(M^,f) relation 
( |Renzini| [20091 |Peng et al.||2010| jBouche et aLl[20T0 
MA 10). Our zC-SINF sample includes by selection 
only star- forming galaxies at 1.4 < z < 2.5, with very 
young stars boosting the Ha emission line. We then 
adopted two different types of SFH: (i) CSFR with an 
age lower limi t of 0. 1 Gyr, usin g stel lar population mod- 
els from b oth ]Maraston] ( |2005| ), and Bruzual and Char- 
Tot] ([2003, hereafter MA05, and BC03, respectively), 
and (ii) SFR exponentially increasing with time (SFR 
oc e +t / r ), i.e., the so-called "inverted r-models" (inv-r, 
cf. MA10), with fixed forma tion redshift (z/ — 5, i.e. 
fixed ages, see Section |3.2.3| ), then using MA05 stel- 
lar population models. As MA05 templates (bo th CSFR 
and in v-r) are built with the lKroupa ( [200 1| ), or Salpeter 
( |1955| ) IMF, we used the following relations t o correct 
galaxy stellar masses, according to the adopted Chabrier 
p003) IMF: log(M^chabrier) = log(M*,Krou P a)-0.04, and 
log(M* ; chabrier) = log(M^ ;S ai P eter) - 0.23. The same correc- 
tions apply to log(SFR). 



In Appendix A, Figure 20 (entirely published as online- 
only material) shows the oest-fit SEDs that we have ob- 
tained with the three mentioned combinations of mod- 
els and SFHs. The corresponding best-fit parameters are 
listed in Tableffl It is worth noticing that the best-fit SEDs 
obtained with aifferent templates are very similar; those 
obtained with BC03+CSFR and MA05+CSFR models al- 
most fully overlap each other. In the next sub- sections the 
SED best-fit parameters (in particular M* and SFR) ob- 
tained in the three cases are compared to each other, thus 
showing what kind of systematics may affect the quanti- 
ties derived from SED fitting procedures applied to our 
z ~ 2 galaxies. 

3.2.2. Constant SFR models with different population 
synthesis models 

The main difference between MA05+CSFR and 
BC03+CSFR models is that the BC03 population synthe- 
sis models do not include the contribution of the stars in 
the thermally pulsing asymptotic giant branch (TP-AGB) 
phase, which instead in MA05 models contribute a signif- 
icant fraction of the luminosity (especially in the near-IR) 
for stellar population ages in the range between ~ 0.2 and 
~ 2 Gyr. As a consequence, to compensate for the lack 
of this contribution, BC03+CSFR models tend to over- 
estimate galaxy stellar masses and ages with respect to 
MA05+CSFR models for galaxies in whic h stellar popu- 
lation in this age range d ominate (MA05, Maraston et al. 
2006[|Wuyts et al.|2007[ etc.). Moreover, for ages in ex- 
cess of ' ~ 0. 1 Gyr MA05 models tend to have redder B— V 
colors with respect to the BC03 models, which is due 
to differences in the adopted stellar evolutionary tracks 
that affect also other evolutionary phases (e.g., different 



treatment of overshooting in the stellar interiors; see the 
top pa nel of MA05 Figu re 27, and also the top-left panel 
of the Bruzual A 2007 Figure 3). The redder colors of 
MA05 model spectra at given stellar age require younger 
best-fit ages in order to reproduce the observed colors, 
and this drives the results towards lower stellar masses 
and higher SFRs. 

This is illustrated in the left panels of Figure [2] and 
Figure [3] comparing masses and SFRs obtainea with 
the BC03+CSFR and MA05+CSFR models for the zC- 
SINF sample. Note that the ages derived with both 
MA05+CSFR and BC03+CSFR models are indeed very 
low for the majority of the zC-SINF galaxies, and tend 
to cluster to the artificially imposed minimum age (i.e., 
0.1 Gyr, see the previous subsection and Table [41. These 
galaxies are shown as red filled circles in the left panels 
of Figure[2]and Figure(3] and for them MA05+CSFR and 
BC03+CSFR models give consistent results. Instead, re- 
sults differ for galaxies for which at least BC03+CSFR 
ages exceed 0.2 Gyr: MA05+CSFR masses are smaller, 
and SFRs are higher, both by a factor of ~ 1.5, compared 
to those obtained with BC03+CSFR models. These dis- 
crepancies can be ascribed to a combination of the two 
effects mentioned above, MA05 optical colors being red- 
der than BC03 ones for ages older than ~ 0.1 Gyr, and 
near-IR color being redder for ages older than ~ 0.2 Gyr. 

3.2.3. Constant SFR models vs. models with exponentially 
increasing SFR 

An option for avoiding excessively young ages from 
the best fit procedures is to fix the galaxy formation red- 
shift (zu as the beginning of star formation) to some high 
value. The age then follows from the cosmic time dif- 
ference between the formation redshift Zf and the red- 
shift 

£spec a t which the galaxy is observed. Galaxies in- 
deed abound at z > 2 (and even at z ^> 2), and are the 
likely precursors of those we observe at z ~ 2. This lat- 
ter procedure was adopted in MA 10 in conjunction with 
"inverted-r" models, by fixing Zf — 5, and r > 0.3 Gyr, 
and the same recipe is used here, then comparing the re- 
sults to those of CSFR models. The outcome is fairly in- 
sensitive to the precise value of Zf, because most stars are 
formed in the last r time interval in these inverted-r mod- 
els (MA 10). In practice, since most of the galaxies in our 
sample are clustered at z ~ 2-2.5, the assumption that 
star formation begins at z — 5 implies galaxy ages in the 
range ~ 1.1-2.1 Gyr. Clearly, the two objects at z ~ 1.4 
turn out to be older, with ages ~ 3 Gyr. The values of r 
selected by the best-fit procedure range from ~ 0.345 to 
0.89 Gyr. 

Stellar masses and SFRs derived assuming a constant 
SFR and those obtained from inverted-r models (in both 
cases using MA05 stellar population models) are com- 
pared in the right panels of Figure [2] and Figure [3] Stel- 
lar masses from inverted-r models are on average ~ 1.5 
times higher than those obtained assuming a constant 
SFR, and the SFRs are ~ 2 times lower. Therefore, spe- 
cific SFRs (i.e., SFR/M*) are on average a factor ~ 3 
lower when adopting inverted-r models. As discussed in 
the previous section, and reported in Table 4, best fit ages 
are typically very short with CSFR models. Since for 
such models =AgexSFR (apart from the mass return 
fraction) the best-fit procedure picks a higher SFR (and 
a higher Ay) in order to match the M*-sensitive near- 
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Fig. 2. — Comparison among zC-SINF galaxy stellar masses derived by using different stellar population models, and different 
SFHs. Left panel: MA05+CSFR vs BC03+CSFR. Red filled circles are galaxies with both Age (M A05) and Age (B co3) ^ 0.1 Gyr, i.e. 
the imposed minimum age, and open circles galaxies with (at least) Age(BC03) > 0.1-0.2 Gyr (see the discussion in the text). Right 
panel: MA05+CSFR vs MA05+inv-r. Here the red filled circles indicate galaxies with Age(MA05+csFR) — 0.1 Gyr (i.e., artificially 
set to the imposed minimum value), and the open circles galaxies with Age(MA05) > 0.1-0.2 Gyr. It is evident that the best-fit stellar 
masses of the latter better agree with masses estimated using inverted-r models. 
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Fig. 3. — Same as Figure [2] but with the SFR, instead of mass. A comparison with Figure [2] indicates that inv-r models lead to a 
factor no 3 reduction in the specific SFR, compared to CSFR models. 



IR flux. In doing so, the "best" compromise is one in 
which is underestimated and SFR and Ay are over- 
estimated, compared to the inverted-r models, which in- 
stead are allowed to accumulate mass for typically a ~ 10 
times longer time. 

Incidentally, it is worth noting that such a large differ- 
ence in specific SFRs would have an impact on the com- 
parison of the integral of the cosmic SFR d ensity with the 
growth of cosmic stellar mass density (cf. Wilkins et al. 



2008). Note that these systematic differences are much 
smaller for the few objects for which CSFR+MA05 mod- 
els give ages older than 0.2 Gyr (open circles in FiguresE] 
and [3]) compared to those for which CSFR+MA05 mo(> 

els give ages below 0.2 Gyr (red filled circles). 

For the reason mentioned above and in Sections |3.2.1| 
we favor inverted-r models for the description of our 
z ~ 2 actively star-forming galaxies. However, in the fol- 
lowing we keep reporting all three sets of SED-fitting 
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results (CSFR with BC03 and MA05, and inv-r with 
MA05) in order to better gauge the robustness (or the sys- 
tematic uncertainties) of the results. We note in fact that 
the x 2 procedure often gives formal error in M* or SFR 
which are just a few per cent. Such exceedingly small 
formal errors clearly contrast with the much larger sys- 
tematic differences noted above. Given the benchmark 
nature of these SINFONI galaxies, quoting the results for 
various assumed SFHs will also favor intercomparisons 
with future datasets. 

3.3. Comparison with other samples and selection-bias 

In this section we compare properties such as colors, 
masses, and SFRs of the zC-SINF galaxies with those of 
the full sample of K < 23 galaxies at z ~ 2 in COSMOS 
(E. Daddi and H.J. McCracken private communication), 
and with those of the SINS Ha sample of FS09 (i.e., a 
similar sample at the same redshifts). These comparisons 
are merely meant to establish which are the main selec- 
tion biases affecting the zC-SINF sample, and to which 
extent it is representative of the general population of 
galaxies at z ~ 2. Besides the bias relative to the red- 
shift distribution, discussed at the end of Section[2j one of 
the most significant biases affecting the zC-SINF sample 
is due to the zCOSMOS-Deep B AB < 24.75-25.25 cut- 
off, applied to reach an adequate S/N ratio with VIMOS 
spectroscopy (see Section [2j). Besides excluding virtu- 
ally all passively evolving (quenched) galaxies at z > 1 .4, 
this selection criterion also excludes the reddest K < 23 
star-forming galaxies at 1.4 < z < 2.5, i.e., objects highly 
obscured by dust, which are the majority of activel y star 
forming and massive galaxies in this redshift range (Ren- 
zini and Dadd i|2009| ). This is shown in Figure [4] where 
the positions of the zC-SINF galaxies in the BzK dia- 
gram (red filled circles) are compared to those of the gen- 
eral population of galaxi es in the COSMOS field from 
McCracken et al.| ( |2010| ), down to the (70% complete- 
ness) limit K AB < 23. Above the diagonal line, paral- 
lel to the reddening vector, lie the sBzKs, the upper-right 
wedge-shaped region is occupied by the pBzKs, while the 
lower right part of the plot is mostly populated by galax- 
ies at z < 1 .4 and by the stel lar sequence identified by 
the red line ( |Daddi et al.|2004| ). All the zC-SINF sources 
are in the sBzK region, but one, i.e., object ZC403103 
which was selected for VIMOS spectroscopy following 
the UGRfBX criterion (w* -g + ~ 0.74 and g + -r + ~ 0.36). 
Although this object lies in the z < 1.4 galaxy region, it 
has a solid spectroscopic redshift ofz = 2.3613 (Lilly et 
al. 2011, in prep.), and its location in the BzK plot is 
as expected for very young, almost unreddened popula- 
tions at z > 1.4 (Daddi et aLl|2004| ). To better illustrate 
the effect of the B < 25 cutoff used in zCOSMOS-Deep, 
the COSMOS sample is split in two sub-samples, with 
B AB < 25 (black dots), and B AB > 25 (green dots), re- 
spectively. Clearly the zC-SINF galaxies are among the 
bluest sBzKs, with (z~K) and (B-z) colors < 2-2.5, 
whereas redder sBzK (green dots in the top part of the 
diagram) are missed by the zCOSMOS-Deep selection. 
These are the most dust-obscured galaxies, including the 
most massive/star-forming o bjects, which are excluded 
by the B < 25 constraint (see Renzini and Daddi 2009 ). 
We also compare the colors of the zC-SINF sample with 
those of the FS09 SINS objects with measured BzK col- 
ors (43 out of 62), represented by open squares in the fig- 



ure. The SINS galaxies have BzK colors very similar to 
the ones of our sample, because comparable optical mag- 
nitude cutoffs were also applied for the parent spectro- 
scopic samples from which the SINS targets were drawn 
(cf. FS09). FS09 have argued that such optical selection 
excludes more than 5 0% of the z ~ 2 galaxies in the same 
mass range (see also [va n Dokku nTet al.||2006{ |Rudnick| 
|et al.|2006{ |Grazian et al.|2007| ), a loss traction that may 
exceed ~ 80% for gala xies more massive than ~ 1O 1O M 
in the COSMOS field ( [Renzini and Daddi|2009| ). 




2 4 

(B-z) 



Fig. 4.— BzK diagram for the total 2 deg 2 COSMO S sample of 
^-selected s ources (K < 23, McCracken et al. 2010 ). Following 
Daddi et al. (2004), the magenta lines mark off the three differ- 
ent regions of sBzK, pBzK, and low redshift galaxies. The red 
curve in the bottom p art of the diagram sh ows the expected col- 
ors for the stars (from Lejeune et al. 1997), overlapping with the 
locus of real star sequences. The zC-SINF targets are shown as 
red filled circles, and are all included in the sBzK left part of the 
diagram, but one, i.e., an U GR/BX- selected star-forming galaxy 
at z = 2.361, located in the low-redshift area (see discussion in 
the te xt). Black and green d ots are K < 23 COSMOS BzKs 
(from McCrack en et al.|2010) with B AB < 25, and B AB > 25, re- 
spectively. The SINS sample of F09 (open white square) is also 
shown for comparison. 

Another bias is illustrated in the left and middle pan- 
els of Figure [5] where the SFR-M* relation of the zC- 
SINF sample is compared with that of the whole sam- 
ple of COSMOS sBzK galaxies (left panel), and with the 
SINS sample (open squares in the middle panel). While 
our zC-SINF sample covers, by selection, the same mass 
range of COSMOS sBzKs, at the low-mass end of the 
SFR-M* relation, it includes preferentially galaxies with 
the highest SFR due the lower limit imposed to the SFR 
(> 10 M /yr, i.e., see item 4 in Section [2]). On the other 
hand, at the high-mass end of the relation we are missing 
the most star-forming (dusty) objects with B > 25 (i.e., 
the green dots in the top-right part of Fig. [5j left panel). 
The combined effect of these selection criteria flattens the 
SFR-M* relation of zC-SINF galaxies compared to that 
of the whole population at 1.4 < z < 2.5. We note that, 
in general this effect tends to appear whenever a selec- 
tion criterion sensitive to the current SFR is applied on 
a sample of star-forming galaxies previously selected by 
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mass. In fact, especially at low values, galaxies with 
SFR above the average SFR-M* relation will be prefer- 
entially selected, a typical Malmquist bias. The middle 
panel of the same figure shows that the same effect is 
seen also in the SINS Ha sample (open squares; see also 
FS09), that was selected with similar criteria, and then 
suffers the same bias of our zC-SINF sample. All the 
SFRs reported in Figure [5] are derived from the UV rest- 
frame luminosity ( SFR^y), according to the relation of 
Dadd i et al.| ([2004), and ad justed to a Chabrier IMF (as 
explained in Section |5.1.1| ). The SFRf/y derived for the 



-0.006, and la dispersion a 



log 



M* K 



M*!- C AT 



= 0.127, and 



zC-SINF sample are also listed in Table [6] The stellar 
masses shown on the left panel were determined using the 
e mpirical relation bas ed on the K magnitude, also derived 
by paddi et al.| ( [2004] ), so to ensure a homogeneous com- 
parison between zC-SINF galaxies and COSMOS sBzKs, 
for which masses from SED fitting were not available for 
the full sample. As the paddi et aL| ( |2004| ) relation was 
calibrated on stellar masses derived from SED fitting with 
BC03 models and Salpeter IMF, we applied the appro- 
priate corrections to Chabrier IMF. On the contrary, all 
the stellar masses reported in the middle panel are from 
SED fitting with BC03+CSFR models, to avoid systemat- 
ics derived from different assumptions in the comparison 
with the SINS sample (open squares). The right panel 
of the same figure shows the comparison between masses 
computed with the two methods, which are on average in 
good agreement, with a scatter of ~ 0.2 dex. 

As a result, the zC-SINF sample offers an adequate set 
of benchmark galaxies for 2 < z < 2.5, with a few se- 
lected objects probing later epochs (z ~ 1.5). 

3.4. Comparing results from different catalogs 

A final check is relative to the possible impact on de- 
rived quantities (i.e., M+, and SFR) of the use of differ- 
ent photometric catalogs for galaxies in the COSMOS 
field. We took advantage of the public COSM OS /-band 
selecte d catalog (hereafter, / - CAT) from jllbert et aL 
2009) to compare the results obtained using the -band 



selected catalog (K-CAT) ofjMcCracken et al.|(2010), 



and therefore quantify the differences in physical quan 
tities when using these two catalogs. The I -CAT dif- 
fers from the K-CAT adopted in this work mostly in 
: (i) the optical selection (I AB < 26), and source extrac- 
tion at the /-band positions; (ii) the PSF-matching: all 
the images (from w* to ^-band) were degraded to the 
same PSF of l."5 (i.e., the /T-band PSF) as detailed in 
|Capak et aL] ( |2007[ ), while no PSF degradation was used 
in the K-CAT; (ni) the magnitudes measurement in 3" 
diameter apertures (instead of 2"); (iv) the I -CAT does 
not include //-band photome try. We used the same pro- 
cedure described in Section |3.2| to fit the galaxy SEDs 
with the I -CAT photometry, using the HyperZmass soft- 
ware, with BC03+CSFR models and Chabrier IMF. We 
also computed the SFR from the UV rest-frame luminos- 
ity (i.e., in particular, from the observed g-band magni- 
tudes), applying the s ame re lation used for the K-CAT 
(described in Section [5. 1.1| ). Then, we compare the re- 
sulting masses and SFRs with those obtained using the 
K-CAT, and show the results in Figures [6| [7] and [8] 

Stellar masses from SED fitting and SFR from UV 
derived using such different catalogs are generally in 
good agreement with each other, with an average off- 

set of /j \log^A = 0.009, and /j \log SFR 



a 



log 



SF Ryy^K-CAT 
SFR(JV,I-CAT 



= 0.132, respectively. The correspond- 
ing individual differences and their distributions are 
shown in Figures [6| and [7] In addition to the system- 
atic offset, one notes an apparent trend in the ratio of 
SFR^/y estimates as a function of SFR UViK - cat . We in- 
vestigated the possible cause of this trend, and found that 
it is due to a combination of two effects. First, system- 
atic fainter 5-band magnitudes from the K-CAT with re- 
spect to the g magnitudes from I -CAT (with a mean off- 
set of < B k-cat ~ 8i-cat >~ 0.37 mag) lead to a smaller 
K-CAT UV rest-frame luminosity (Lisoo)- Second, dust 
extinction parameters, A 1500, derived from I -CAT are 
smaller than those from K - CAT for the most obscured 
sources (which are also the most star-forming galaxies), 
while the values are comparable for the less obscured 
sources. 

On the other hand, when using SFRs from SED 
fits, a slight systematic shift is present between 
the SFRs derived from the two different catalogs, 



i.e., fi 



log 



log 



SFRse. 



SFR SEj 



= 0.065, with a dispersion 



SFRi-cat 



SFRjjv,I-CAT 



- 0.079. This is due to systematic smaller 

values of the E(B-V) parameter derived fitting SEDs 
with the I -CAT compared to those obtained from the 
K-CAT. The corresponding individual differences and 
their distribution are shown in Figure [8] 

In summary, these comparisons offer an estimate of the 
reliability of the derived masses and SFRs of zC-SINF 
galaxies here reported in Table [4] By taking into account 
the results of these sections, we can state that zC-SINF 
galaxy masses and SFRs derived from SED fitting can be 
considered reliable within ~ 0.2-0.3 dex (i.e., within a 
factor of 1.5-2.). 

4. SINFONI DATA 

4.1. Observations 

SINFONI ( psenhauer et al.|2003a[|Bonnet et al.|2004] ) 
is a near-infrared integral field spectrometer, mounted 
at the Cassegrain focus of UT4 at the Very Large Tele- 
scope (VLT). The instrument is a combination of the inte- 
gral field spectrograph SPIFFI (Eisenhaue r et al.|20 03b) 
and th e MACAO adaptive optics module ( |Bonnet et al.] 
2003). SPIFFI slices the square field of view (FoV) into 
32 slitlets and rearranges them into a one-dimensional 
pseudo long-slit, which is then fed to the grating spec- 
trometer part of the instrument. Each of the 32 slitlets is 
projected onto 64 detector pixels, thus the spatial pixels 
(called spaxels) on the sky are rectangular. The angu- 
lar size of the slitlets can be selected between 250, 100, 
and 25 mas (corresponding to pixel scales of 125, 50, and 
12.5 mas, and FoV of 8" x 8", 3" x 3", and Of'8 x 0f'8, 
respectively). Spatial dithering of on- source exposures by 
an odd or fractional number of pixels during the observa- 
tions ensures adequate sampling of the spatial axis per- 
pendicular to the slitlets. In this way, SINFONI records 
simultaneously the spectrum of every spaxel across the 
FoV. In the reduction procedure, the three-dimensional 
datacube is then reconstructed, with two spatial dimen- 
sions and one spectral dimension. SINFONI can be used 
with or without the AO module, and in the latter case the 
spatial resolution is dictated by the natural seeing. 
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Fig. 5. — Left and middle panels: The SFRuv-M* relation found for the zC-SINF targets (red filled circles) is compared with that 
of the COSMOS sBzKs (dots in the left panel), and of the FS09 SINS sample (black open squares in the middle panel). In the left 
panel the black dots are the COSMOS sBzKs with Bab < 25, and the green ones are sBzKs with Bab > 25 (as explained in the text). 
The red solid line show the best-fit of the SFRuv-M* relation for COSMOS sBzKs, also reported in the middle panel for comparison. 
The zC-SINF targets have quite high SFRs compared to the total sample of sBzK, especially at low masses, and cover the full mass 
range of K < 23 selected sBzK. The resulting SFR-M* for the zC-SINF sample is therefore appreciably flatter than that relative to the 
full COSMOS sample. Note also the zC-SINF galaxy with the lowest mass (the faintest TT-band magnitude, i.e., ZC403103), which 
was selected according to the U GR criterion. Right panel: Comparison between zC-SINF stellar masses estimate d with two differen t 
methods: (i) log(M*)zx)4, used in the left panel, derived from the empirical relation based on the K magnitude (Daddi et al. 2004 ); 
(ii) log(M*)sED, used in the middle panel, derived through SED fitting (BC03+CSFR models). The solid and dashed lines show the 
best-fit of the relation between these quantities, and the 1 a scatter, respectively. 
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Fig. 6. — Left panel: Comparison between stellar masses de- 
rived for the zC-S INF sample by fitt i ng SE Ds ( BC03+CSFR 
model s) built with |McCracken et al.| ( |2010| >, and |Ilbert et alT] 
(2009 ) photometry (\og{M) K - C at, and \og{M)i- cat , respectively). 
Right panel: the histogram shows the mass -difference distri- 
bution of the sample, together with the corresponding Gaussian 
fit. On the top part of the diagram, the number of galaxies used 
in the fit (#), the mean (/x), and the standard deviation (a) of the 
distribution are labeled. 



The SINFONI data presented in this paper were mostly 
taken as part of an ESO VLT Large Program (PI: A. 
Renzini), which started in February 2009. As already 
mentioned, our sample also includes objects observed in 
earlier SINFONI programs (see Table IT]): four "pilot" 
zCOSMOS targets observed as part ofthe MPE SIN- 
FONI Guaranteed Time Observations (GTO; data taken 
in spring/summer 2007, PI: N.M. Forster Schreiber), as 
well as the five zCOSMOS targets from a Normal Pro- 
gram (PI: A. Renzini) observed in February 2008. The 
targets for all these three programs were selected and ob- 
served in the same fashion, and with identical scientific 
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Fig. 7. — Sa me as Fig.[6[ but with SF Rt/y instead of mass. The 
SFRf/y from |McCracken et al. | pOTO) i.e., log(SF R) K - cat ) was 
derived from the B-band magnitudes, and that from Ilber t et aL] 
( |2009[ i.e., log(SFR)i- ca t) was derived from g-band magnitudes 
as explained in the text. 

goals. They are therefore presented and discussed to- 
gether in this and future papers of this series. In this 
paper, we focus on the results for the first part of the 
program and present the analysis of the no- AO data for 
the full zC-SINF sample. While for the AO observations 
the intermediate SINFONI pixel scale, i.e. 50 mas, will 
be used to achieve the full gain in resolution afforded 
by AO, for the no-AO data the largest SINFONI pixel 
scale of 125 mas has been chosen. The reason is pri- 
marily for the higher surface brightness sensitivity, and 
for a better optimization of the AO follow-up observ- 
ing strategy. Such optimization strategy includes point- 
ing, field rotation, strategy for the sampling of the back- 
ground emission (appropriate with the smaller FoV used 
in AO mode), and proper estimate of the total integration 
time required to reach sufficient S/N per resolution ele- 
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Fig. 8. — Same as Fig.[6]and Fig. 17] but with SFR derived from 
SED fitting (BC03+CSFR models). 

ment to achieve our AO science goals. The AO integra- 
tion time is based on the no- AO Ha surface brightness 
distribution of each individual galaxy, and ranges from 
about 4 hours to 10 hours. We refer to the forthcom- 
ing AO papers for the exact details for each target. We 
used the H or K band gratings in order to map the rest- 
frame 6562.8 A Ha emission line for galaxies in the red- 
shift ranges 1.4 < z < 1.7 and 2 < z < 2.5, respectively 
(also taking into account the OH avoidance). With these 
choices of grating and pixel scale, the nominal spectral 
resolution is R ~ 2900 and 4500 in the H and K bands, 
respectively. 

The natural- seeing observations were carried out in 
Service Mode between February 2008 and April 2010 
(see TablefTT). For each object a typical science "observ- 
ing block' ^OB) included six on- source dithered expo- 
sures, each one with exposure time of 600 s, for a total 
OB integration time of 1 hr. For each pair of successive 
exposures, the target was positioned in opposite quad- 
rants or halves of the FoV; additional small dither offsets 
of « 1/10 of the full FoV were applied to avoid redun- 
dant source positions on the detector frame. With this 
strategy, the area with full integration time (hereafter the 
"effective FoV") ranges from « 4" x 4" to w 4" x 8", 
depending on the source. Exposures of the acquisition 
stars used for blind offsetting to the galaxies were taken 
just before the science target data to monitor the seeing 
and the positional accuracy. Standard stars for absolute 
flux calibration and correction for atmospheric transmis- 
sion (late-B, early- A, and G1V to G3V stars, with near- 
IR magnitudes of 7-10 mag) were observed with the same 
instrument setup, at similar airmass, and close in time to 
the science targets. 

4.2. Data reduction 

We reduced the data with the SPRED softwa re, specif- 
ically developed for SINFONI observations (Schreiber 
et al.||2004[ |Abuter et aL| |2006 ), along with additional 
custom routines. We followed the procedure optimized 
for data of faint, high-redshift targets described by FS09, 
where full details can be found. Notably, this pr ocedure 
features the method developed by Davies (200>7]) for ac- 
curate wavelength calibration and background subtrac- 
tion, which substantially reduces residuals from night sky 
emission lines. Each individual exposure is reduced sepa- 
rately and reconstructed into a wavelength-calibrated and 



distortion-corrected three-dimensional data cube. The fi- 
nal reduced cube of each science target is obtained by 
combining the spatially-aligned individually reduced data 
cubes. We combined the cubes by averaging them with a 
sigma-clipping algorithm (typically clipping at the 2.5 a 
level). This step also produces a "sigma-cube," whose 
pixel values correspond to the standard deviation of pix- 
els at the same spatial and spectral coordinate in the cubes 
for individual exposures (see Section 4.2 and Appendix C 
ofFS09). 

The data of the acquisition and telluric standard stars 
were reduced following the same procedure as for the 
science data. For each target, the spatial point-spread 
function (PSF) was measured on the two-dimensional im- 
age of the acquisition star, obtained by averaging (with 
a-clipping) all the wavelength channels of the star's re- 
duced data cube. The stellar images were fit with Gaus- 
sian profiles, which represent adequately the SINFONI 
spatial PSF in no-AO mode (see Appendix B of FS09). 
The effective angular resolution of our data sets from the 
PSF full width at half-maximum (FWHM) ranges from 
0"29 to 0."84, with a mean and median of 0"60 (corre- 
sponding to ^2.5-7 kpc and « 5 kpc, respectively, at 
Z ~ 2), and is listed in Table [T] together with informa- 
tion about the observing periodTthe used grating, and the 
total integration time (T^[s]). The conversion from in- 
strumental counts to flux units were derived by compar- 
ing the integrated star counts from the synthesized broad- 
band image of the star to its 2MASS catalog magnitude. 

The effective spectral resolution of the reduced data 
was measured from unblended, bright night sky lines 
based on "sky" data cubes created in a similar way as 
the reduced science data cubes but without background 
subtraction. The night sky line spectra were extracted by 
integrating over a ~30 pixel-wide square aperture. In re- 
duced, seeing-limited data taken with the 125 mas pixel 
scale, the night sky line profiles are well approximated by 
a Gaussian function and show little dependence on spa- 
tial position. The FWHM in wavelength units is constant 
across each band and, for our data, corresponds to « 120 
and « 85km s -1 in H- and K-bznd, respectively. 

4.3. Measurements of the Ha emission line properties 
4.3.1. Integrated Ha properties 

We detected the Ha line emission for 25 of the 30 
galaxies observed. For each detected target, we mea- 
sured the global Ha flux, systemic redshift, and veloc- 
ity dispersion from spatially-integrated spectra in circu- 
lar apertures centered on the centroid of the line emission 
determi ned fro m the line maps (extracted as described in 
Section [4. .3. 2\ . The aperture size was chosen to enclose 
the total flux estimated from a curve-of-growth analysis 
within the deepest part of the effective FoV. Typical aper- 
ture radii (r ap ) for our objects range from 0."75 to 1"5, 
and are reported in Table [5] 

The Ha properties were derived using the code LINE- 
FIT (devdope^yJl^INS team and described by FS09 
and |Davies et al.||2011| ) assuming that a single Gaussian 
profile represents the intrinsic emission line profile of 
the source. In the fitting procedure, LINEFIT takes into 
account the instrumental spectral resolution by convolv- 
ing a Gaussian profile with the spectral response func- 
tion rep rese nted by empirical night sky line profiles (see 
Section |4~2| ). The resulting profile is fit to the observed 
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spectrum by means of a minimization procedure that ac- 
counts for the noise spectrum of the data through Gaus- 
sian weighting, and the formal uncertainties of the best- 
fit parameters are derived from 100 Monte Carlo simula- 
tions. The noise spectrum is computed fro m th e "sigma 
cube" associated with each data set (Section |4~2] ) and fur- 
ther accounts for deviations from Gaussian scaling with 
aperture size (due to, e.g., correlated noise) from an anal- 
ysis of the effective noise properties of each data cube 
(see Section 5.2 and Appendix C of FS09 for details). A 
linear continuum component underneath Ha was simul- 
taneously fit with the emission line profile; in most cases, 
however, the continuum is essentially undetected in our 
~ 1 hr integration SINFONI data. 

The total Ha fluxes F(Ha) are given in Table |5j along 
with the vacuum redshifts z(Ha) and integratea veloc- 
ity dispersions cr(Ha)i ntegrated corresponding to the cen- 
tral wavelength and width of the best-fit Gaussian profile. 
For undetected sources, the 3 a upper limits on the total 
Ha fluxes are reported. These upper limits were derived 
based on the noise spectrum for a fixed circular aper- 
ture of \" radius centered at the expected position of the 
source, assuming that the Ha redshift is equal to the red- 
shift from the zCOSMOS optical spectroscopy and that 
the intrinsic line width is cr(Ha)i nte grated = 130 kms~ l (the 
average value for the detected sources). We emphasize 
that the velocity dispersions listed in Table [5] and quoted 
throughout this paper are corrected for the instrumental 
velocity broadening. We also note that, given the mea- 
surement procedure, the global cr(Ha)i ntegrateci values in 
Table [5] include a contribution from internal velocity gra- 
dients (e.g., due to rotation), as they are obtained from the 
integrated spectra without any shift of individual pixels to 
the systemic velocity. The uncertainties from the absolute 
flux calibration are estimated to be ~ 10% and those from 
the wavelength calibration, < 5%. Together with uncer- 
tainties from the continuum placement an d the ch oice of 
aperture for the measurements (see Section 4.3.31) , we es- 
timate that systematic uncertainties amount to 30% typi- 
cally, and up to 50% for the objects with the lowest S/N 
ratio. 

Table[5]also lists the fractional contribution frac B B(Ha) 
from the integrated Ha line flux to the total broadband 
flux density from the observed magnitude in K or H band 
(Table [3} for the targets at z > 2 or z < 2, respectively. 
The associated uncertainties reported in the Table corre- 
spond to the formal measurement uncertainties of the Ha 
flux and broadband magnitudes. The combined system- 
atic uncertainties (including those of the absolute flux cal- 
ibration) are largely dominated by those of the SINFONI 
Ha data and are thus typically 30%, and up to 50% for 
the faintest objects. 

4.3.2. Ha emission and kinematic maps, position-velocity 
diagrams, and sizes 

We extracted maps of the velocity-integrated Ha line 
flux, velocity (relative to the systemic velocity), and ve- 
locity dispersions for the zC-SINF sample in the same 
manner as described above for the integrated properties, 
but now applying the procedure to the spectrum of in- 
dividual spaxels. Prior to the fitting, the data cubes 
were lightly smoothed with 3 pixel-wide median filter- 
ing in both spatial and in the spectral dimensions, in 
order to improve the S/N ratio. Because of the short 



« lhr integration time for our pre-imaging SINFONI 
data, the kinematic maps for a large fraction of the sample 
are significantly affected by low S/N even after median- 
filtering, and hamper reliable extraction of the full two- 
dimensional information. Higher S/N as well as higher 
spatial resolution will be achieved for the subset of the 
sample followed-up with SINFONI+AO observations. 

For the pre-imaging data sets, useful kinematic in- 
formation can nevertheless be obtained from position- 
velocity (p-v) diagrams, as well as profiles of the ve- 
locities and velocity dispersions extracted along the ma- 
jor axis of the objects. For the p-v diagrams, we used 
a synthetic slit with width of 6 pixels (0"75) in all cases, 
integrating the spectra along the direction perpendicular 
to the slit orientation. The slit was always positioned so 
as to pass through the centroid of the Ha emission based 
on the extracted line map. Whenever possible, the slit 
was aligned with the kinematic major axis identified from 
the velocity fields, and defined as the direction of steep- 
est observed velocity gradient across the source. In those 
cases, the Ha kinematic major axis is generally consis- 
tent with the morphological major axis identified from 
the Ha line maps as well as the HST ACS /-band imag- 
ing. In the other cases (i.e., for data with too low S/N), 
the slit was aligned with the morphological major axis, 
and we verified with p-v diagrams extracted along sev- 
eral different P.A.s that larger velocity differences were 
not measured for different directions than along the mor- 
phological P.A. The maximum observed velocity differ- 
ence, V bs = V max — 

Vmin, was derived from the major axis 
velocity profile obtained from fits to spectra integrated in 
contiguous, 6 pixel- wide apertures equally spaced along 
the major axis (with centers separated by 3 pixels). Sig- 
nificant velocity differences could be measured in 17 of 
the 25 detected sources from the pre-imaging no- AO data 
sets presented in this paper. 

To estimate the Ha sizes of the detected zC-SINF ob- 
jects, we determined the half-light radius based on the 
curve-of-growth analysis in circular apertures described 
in Section |4.3| The intrinsic sizes r^Ha) are cor- 
rected for the effective spatial resolution by subtract- 
ing in quadrature the PSF half width at half-maximum 
(HWHM) appropriate for each data set. For one source 
(ZC4 15087), the observed half-light radius is smaller 
than the PSF HWHM; since the PSF was not observed 
simultaneously with the data, this could be attributed to 
variations in time of the seeing between the observations 
of the acquisition star used as PSF calibrator and the sci- 
ence target. For this galaxy, we adopted the observed size 
as upper limit to the intrinsic size. The formal uncer- 
tainties of the r!/ 2 (Ha) were estimated as in FS09 (and 
turn out to be very similar as well), accounting for vari- 
ations of ~ 20% in the PSF FWHM between successive 
OBs and for deviations from axisymmetry of the PSF pro- 
files (the median ellipticity of the PSFs associated with 
our data sets is ^0.1); the resulting uncertainties on the 
r!/ 2 (Ha) measurements are typically w 35% (median). 
These uncertainties do not account for other potentially 
important factors affecting size determinations, such as 
the dependence on the sensitivity of the data, the actual 
intrinsic surface distribution of the Ha emission, and the 
measurement method. We address these issues further in 
the next subsection. 

The velocity-integrated Ha line maps and the position- 



The SINS-zCOSMOS Project 



13 



velocity diagrams of all detected objects in the zC-SINF 
sample are shown in Appendix B, Figure [21] (entirely 
published as online-only material), along witn the ACS 
/-band images and the source-integrated spectra around 
Ha. The maximum observed velocity difference mea- 
surements are given in Table [5] 

4.3.3. Impact of the surface brightness sensitivity on the 
measurements of the Ha properties 

In view of the relatively short integration times of about 
1 hr for the zC-SINF pre-imaging SINFONI data sets, we 
examined the possible impact of surface brightness sen- 
sitivity on the derived total Ha fluxes and half-light radii. 
For the measurements presented in this work, we chose 
a curve-of-growth analysis in circular apertures for its 
simplicity, and because it does not rely on assumptions 
about the intrinsic spatial distribution of the line emis- 
sion. This approach also allowed us to derive in a uniform 
and consistent manner size estimates for all detected ob- 
jects, which span a significant range in apparent sizes and 
were observed under different seeing conditions. 

We point out, however, that other methods would lead 
to more robust size estimates for the better resolved 
sources with higher S/N data, as discussed by Bouche et 
al. (2011, in preparation; see also FS09). For instance, 
the use of elliptical apertures to construct the curves-of- 
growth would be more appropriate for sources with ellip- 
tical isophotes, such as inclined disk-like systems. Our 
use of circular apertures gives r^Ha) values that are 
analogous to circularized half-light radii; the latter are by 
definition smaller than the actual half-light radius by a 
factor of y/b/a, where b/a is the ratio of the projected 
size along the minor and major axes of the galaxy. In 
addition, and more importantly, the curves-of-growth are 
affected by the surface brightness sensitivity of the data. 
Estimating the impact on the total flux measurements is 
complicated by the a priori unknown details and extent of 
the intrinsic spatial distribution of the Ha line emission. 

To gauge the impact of surface brightness sensitivity on 
the total line fluxes, we used data of the SINS Ha sam- 
ple with a range of integration times from typically a few 
hours and up to lOhr per object. The curve-of-growth 
analysis for the 40 SINS objects with integrations longer 
than 1 hr (with mean and median of 4 hr) was repeated 
after adjusting the surface brightness limits of each data 
set to the equivalent limits for 1 hr integration time. More 
specifically, we determined the average S/N per pixel at 
the radius r ap at which the curve-of-growth converges, 
based on the data with full integration time 7] nt , and used 
to compute the total Ha flux. To mimic the shallower sur- 
face brightness limits for 1 hour integrations, we scaled 
the above average S/N per pixel by 1 / y/T^ appropriate 
in the background-limited regime and identified the new 
r ap for 1 hour as the radius at which the S/N per pixel cor- 
responds to the scaled value. From the curve-of-growth 
and total flux measured within this new radius, we re- 
calculated the ri/ 2 (Ha). This is clearly a simplistic ap- 
proach, but we verified for a few objects that consistent 
results would be obtained if we had used cubes combin- 
ing only a subset of the exposures amounting to 1 hour in- 
tegration time to estimate the r!/ 2 (Ha) in shallower data. 

The results indicate that the total flux derived from the 
shallower data is on average ~ 10% lower than for a typ- 
ical integration time of 4 hr, albeit with significant scat- 



ter. Unsurprisingly, the median differences tend to in- 
crease for data sets with higher integrations, up to ~ 30% 
for 7- 10 hr. The scatter also increases, probably as a 
result of the diverse source morphologies. This exer- 
cise further indicates that the half-light radii could be un- 
derestimated by comparable amounts of ~ 10% -30%. 
This simple analysis thus suggests that while undoubt- 
edly present, surface brightness limitations in our zC- 
SINF pre-imaging data do not appear to lead to substan- 
tial (i.e., by factors of two or more) underestimates of the 
total Ha fluxes or sizes of the targets, at the very least 
in a comparative sense with the deeper SINS Ha sample. 
The longer integrations of the AO-assisted SINFONI ob- 
servations planned for the second part of our program is 
expected to provide better constraints, in particular on the 
Ha sizes. 

5. CONSTRAINTS FROM Ha PROPERTIES 

5.1. Ha Luminosity, Dust Extinction, and Star 
Formation Rate 

The observed Ha luminosity L ^(Ha) of each galaxy 
was then derived from the Ha fluxes extracted from the 
spatially integrated spectra, neglecting any B ai mer ab- 
sorption from the underlying stellar population (Brinch- 
mann et al.|2004| FS09). HII regions suffer dust extinc- 
tion, hence to find the intrinsic luminosities we correct 
L fo(Ha) for the radiation absorbed by dust. In general, 
the most direct estimate of the amount of dust extinc- 
tion in the nebular regions is based on the recombina- 
tion lines ratio Ha/H/3. In this work we did not observe 
the H/3 emission line, and therefore we used as alter- 
native estimate the extinction based on the color excess 
(i.e., E(B-V)) derived from the galaxy SED modeling, 
as well as from the observed B-z color, as described 
in wh at follows. By studying local galaxies, several au- 
thors dFanelli et al |fT988l |Mas-Hesse and Kunth|[T999l 
May yaeT~al.||2004[ [Cid Fernandes et al.||2005[ |Calzetti 
|et al.|1994 , 2000 ) found that the B aimer emission lines 
H II regions are more extincted than the UV stellar emis 
sion (estimated through SED fitting) by an amount em 
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pirical ly given by the relation proposed by Calzetti et al. 
( [T994l|2000| ): 
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In agreement with this prescription, FS09 found that 
nebular regions are about two times more extincted than 
the stellar continuum in their SINS sample of star- 
forming galaxies atl.3<z<2.6. 

The first direct measurements of the B aimer decrement 
(Ha/H/3 flux ratio) in small z ~ 1.5-2.5 samples provide 
further support for extra attenuation towards the HII re- 
gions relative to the bulk of the stellar population that 



dominates the broad-band cont inuum light ( Yoshikawa 
et al.|2010[[Muzzin et al.|2010[ Buschkamp et a l. 2011, 
in prep.). On the other hand,|Reddy et aL](|2010|) derived 



similar extinction for Ha fluxes and UV continuum for 
UV-selected objects at the same redshifts (see also Sec- 
l 5.1.1| ). Thus, in this work we consider both possibili- 



tion 



ties to correct for dust extinction, and, following the nota- 
tion of FS09, we denote as L°(Ha) and L 00 (Ha) the lumi- 
nosity corrected with A v = A v ,sed, and Ay = A v ,hii (from 
Eq. 0), respectively (i.e., L(Ha) = L obs (Ha) x 10 33Av , 
see Table [6). 
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log (SFR [M /yr]) SED log (SFR [M /yr]) SED 

Fig. 9.— The Ha luminosities (in erg s -1 ) L(Ha) vs the SFR obtained from SED fitting with BC03+CSFR models) for our zC-SINF 
sample (red filled circles) and SINS sample (black open squares). Ha luminosities of sources undetected in SINFONI are plotted at 
their 3 a limits (cf. Table[6|. Left panel: Ha luminosity corrected for dust extinction Ay identical to that derived from the SED fitting. 
Right-panel: the same relation with enhanced extinction correction as indicated. In both the panels the black dashed line shows the 
SFR-Ha luminosity relation from Brinchmann et al. (2004) for a Chabrier IMF. The best-fit and zbl a scatter of our data-points are 
shown by the red solid line, and the red dashed lines, respectively. It is evident that in the right panel the best-fit slope better agrees 
with the one-to-one relation derived from Equation [2] 
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Fig. 10. — The same as in Figure[9] but using MA05+CSFR stellar population models. Here we show the results just for the zC-SINF 
sample, as the SFRs from SED fit are only given for BC03 models for the SINS FS09 Ha sample. 
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The SFR [M© yr" 1 ] is direc tly related to th e Ha lu- 
minosity [erg s -1 ] through the |Kennicutt ( 1998 ) relation, 
that is: 



log(SFR(Ha)) = log(L(Ha))-41. 10-0.23, (2) 

once the correction for the different IMF used in this pa- 
per (i.e. , -0. 23 dex, from Salpeter to Chabrier IMF; see 
Section |3.2| ) has been applied. The resulting Ha lumi- 
nosities and SFRs are then reported in Table [6| 

5.1.1. Comparison among different Star Formation Rate 
indicators 

One of the main advantages of our zC-SINF sample 
is the availability of both Ha integral field spectroscopy, 
and multi- wavelength photometric information. This al- 
lows us to compare the derived Ha properties with those 
estimated from the broadband photometry. In particular, 
in the following sections we compare the SFRs from the 
Ha intrinsic luminosities with those obtained from the 
SED fits and from the UV rest-frame luminosity. 

SFR FROM Ha VS SFR FROM SED FITTING 



In Figures [9)[To| and[TT]we plot L°(Ha), and L 00 (Ha), 
derived in the previous section, as a function of the SFR 
estimated by fitting the galaxy SEDs with BC03+CSFR, 
MA05+CSFR, and MA05+inv-r templates, respectively 
(SFRsed)- From these figures it is evident that for the 
zC-SINF sample (red filled circles) the slope of the best 
fit SFRsed -L(Ha) relation better agrees with the one-to- 
one relation derived from Equation [2] (black dashed line) 
when Ha luminosities are correctea for dust extinction 
with A v ,hii (Equation [T]) instead of A v ,sed- This state- 
ment remains true also taking into account the relatively 
small fraction of Ha flux missed by these 1 hr integr ation 
time SINFONI data (i.e., - 10-30%), see Section [51] 
Our results show that this is the case regardless of the 
stellar population models adopted or the SFH assumed. 
This result is in agreement with the finding by FS09 for 
the SINS targets (overplotted as open squares in Figure [9 
for compa rison), an d in contrast to the conclusions by [Era 
et al.| ( 2006b ) an d Re d dy et al.|p010| ). As already men- 
tioned m Section[5.1| these otner studies of z ~ 2 galaxies 



claimed that in their sample Ha regions are as extincted 
as the stellar UV continuum, and no further correction is 

n eeded for L(Ha ). 

|Reddy et al.| ( 2010| ) suggested that this discrepancy 
could be due to the lower stellar masses and SFRs of their 
sample with respect to the bulk of the SINS sample, and 
that it is possible that the nebular reddening may on aver- 
age be larger for galaxies that are forming stars at a higher 
rate. This lower reddening for lower mass galaxies may 
be related to these galaxies naturally having lower metal- 
licities. Alternatively, it could also be due to the presence 
of a more important olde r and, on average, le ss attenuated 
stellar population (e.g., |Epinat et~aL]|2009| ). Our results 
are consistent with this explanation, as our zC-SINF sam- 
ple includes only objects comparable to those with the 
highest masses and Ha luminosities (hence SFRs) of the 
SINS sample. From Figure [9] it appears that the SINS 
objects (open squares) with lower SFRs and masses do 
not necessarily need a larger dust extinction correction 
to match the L(Ha)-SFRsED relation (dashed line), as in- 
stead needed for the most massive and active SINS galax- 
ies. 



SFR FROM Ha VS SFR FROM UV 

In star- forming galaxies, the rest- frame UV luminos- 
ity over the spectral range of 1500-2800 A is dominated 
by the contribution of young (OB) stars, and does not 
suffer much contamination by the older stellar popula- 
tion. Hence, the h u luminosity over this wavelength 
range, corrected for dust extinction, scales linearly with 
the SFR and can be used as an efficient SFR indicator. 
Based on synthetic stellar population models, different 
calibrations have been obtained to convert UV luminos- 
ity to SFR, depending on the reference wa velength used, 
templates, and IMF (see Kennicutt 1998, and reference 
therein). In this work, we used the relation relat ive to 
1500 A, an d adjusted for BC03 models (following [Daddil 
et al.|2004[ eq uation 5), and for the Chabrier IMF (as m 
Section |3.2.l) : 



log(SFR uv ) = log(L 150 o)-27.95-0.23, 



(3) 



where SFR UV is in unit of [M yr 1 ], and Z4500, in 
[erg s _1 Hz _1 ], is the intrinsic luminosity already corrected 
for dust extinction. As for galaxies at 1.4 < z < 2.5 the 
A=1500 A flux is redshifted in the 5-band (which actu- 
ally samples the 1250-1800 A rest-frame), we used the 
5-band magnitudes and the spectroscopic redshift infor- 
mation to derive the observed 1500 A luminosities. The 
intrinsic 1500 A luminosity, entering equation [3 was ob- 
tained by multiplying this observed luminosity by the fac- 
tor 10°- 4Al500 , where A1500 is the extinction parameter at 
A=1500 A. We computed A1500 from the (B-z) colors 
according the relation: 



^1500 : 



:2.5(£-z + 0.1), 



(4) 

derived from |Daddi elaT] ( |2004| ) Equation 4 (i.e., E(B- 
V) uv = 0.25 (B-z + 0.1) anTAisoo = 10E(B-V) UV , cf. 
Daddi et al. 2007 ), which provides an approximate es- 
timate for the extinction of star- forming BzK- selected 
galaxies at z ~ 2. 

In Figure 12 and Figure [13] we compare the extinc- 
tion parameter E(B -V)uv, derived from the (B-z) 
color, and the SFR from the UV rest-frame luminosities 
(both reported in Table [6} with those estimated by fit- 
ting the galaxy SED with MA05+CSFR, BC03+CSFR, 
and MA05+inv-r models. From these figures it appears 
that SED fits and UV-derived values of E(B-V) and SFR 
correlate quite tightly, with some small offset, < 0.05 
in E(B-V) and < 0.2 dex in SFR. We note that in 
Dadd i et aL| ([2004 ) the relations used to derive both the 
E(B - V) parameter, and the SFRf/y were calibrated us- 
ing BC03+CSFR models. This could be the cause of 
the slight offset in the SFRt/y-SFR(MA05+CSFR) and 
SFRf/y-SFR(MA05+inv-r) relations observed in the left 
and right panels of Figure [13] However, the SFRs de- 



rived with different stellar population models and SFHs 
agree within ~ 0.2 dex with the SFRf/y, and notice that 
the trend is opposite for CSFR and inverted-r models, 
which give systematically lower SFR. 

Figure [14] is similar to Figure [9-11 but shows the 
L°(Ha) and L 00 (Ha) as a function of the S¥R UV , derived 
from Equation [3] rather than from the SED fit. As ex- 
pected, this figure confirms that for our zC-SINF sample 
the HII regions appear to be obscured by roughly twice 
as much compared to the stars. It is worth noticing that 
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Fig. 12. — Comparison between the E(B-V) reddening parameter computed from the B — z colors with the Daddi et aT] p004j > 
relation, and the same parameter derived by fitting the galaxy SEDs with MA05+CSFR, BC03+CSFR, and MA05+inv-r models, left, 
middle and right panel, respectively. 
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Fig. 14. — The Ha luminosity L(Ha) vs SFRuv corrected for different amount of dust extinction. Symbols and lines have the same 
meaning as in Figure [9] 
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the L°(Ha) and L 00 (Ha) shown in this figure were ob- 
tained by using the extinction parameters derived from 
the (B-z) colors instead of the ones derived from the 
SEDs. This was done to minimize the systematic effects 
related to different methods to measure the extinction. In 
fact, we verified that by correcting L b s (Ha) with Avsed, 
obtained using BC03+CSFR models, in Figure 14 we 
would obtain a different slope (0.72 instead of IaM) of 
the best-fit L(Ha)-SFR f/ y relation for the zC-SINF ob- 
jects, and a slightly larger scatter. 

5.1.2. Ha Equivalent Width 

For the 25 Ha-detected objects we computed the rest- 
frame Ha equivalent width (EW(Ha)) as the ratio of the 
Ha line integrated flux to the flux density (expressed 
per unit wavelength) of the stellar continuum. The lat- 
ter quantity was derived from //, or K broad-band mag- 
nitudes, for objects at z < 2, and z > 2, respectively, 
then corrected for the Ha line fractional contribution (i.e., 
frac B B(Ha) in Table [5]). Note that for most of the ob- 
jects fracBB(Ha)< 10%, while only for two galaxies it 
reaches 20%-40% (i.e., ZC410041, and ZC41 1737). The 
resulting EW(Ha)s are reported in Table [6] where for 
the Ha undetections we list equi valent width upper lim- 
its. As for fracBB(Ha) (see § |4.3| ), the uncertainties listed 
for EW(Ha) in the Table correspond to the formal mea- 
surements uncertainties while the systematic uncertain- 
ties, dominated by those of the SINFONI data, are typi- 
cally 30% and up to 50% for the faintest sources. While 
the Ha line integrated flux gives an estimate of the cur- 
rent SFR (as only stars with masses > 10 M and age 
< 20 Myr contribute significantly to the ionizing flux), 
EW(Ha) probes the relative importance of the youngest 
most massive stars relative to the bulk of the stellar popu- 
lation, which is responsible for the continuum emission at 
the same Ha wavelength. Therefore EW(Ha) is a stellar 
age indicator, being sensitive to the ratio of current over 
past (average) star formation rate. 

For galaxies with constant SFR over their lifetime, the 
equivalent width should be in inverse proportion to the 
galaxy stellar age, according to a power law, as the stel- 
lar continuum steadily increases with time, while the Ha 
flux remains almost the same (see also van Dokku m et al.| 
|2004| ). In fact, some authors suggested to derive con- 
straints on galaxy star formation histories b y studying 
EW(Ha ) as a f unction of the age (e.g. FS09, Erb et al. 



2006b see also |Kriek et al.|[20TT| ). However, ages esti- 
mated from SED fitting strongly depend on the shape of 
the assumed SFH, and should be considered with some 
caution, also due to th e parameter degeneracy affecting 
the best-fit choice (cf. Bolzonella et al.||2000). In fact, 



as discussed in Section [3T2| for our actively star-forming 
galaxies at z ~ 2, when using constant SFR models, ages 
are underestimated due to the dominant contribution of 
the youngest stellar population to the stellar light, and 
cluster around the imposed age lower limit. On the other 
hand, ages estimated with inverted-r models seem to be 
more realistic for these targets, but they are imposed by 
the choice of the formation redshift, z/. 

Hence, we examined the equivalent width of our zC- 
SINF sample as a function of the (z-H) color, i.e., a 
quantity tightly related to galaxy age, sampling the rest- 
frame 4000 A break for galaxies at z = 1 .4-2.5, and inde- 
pendent of the stellar population models employed (Fig- 



ure [15]). However, as the (z-H) color depends also on 
the dust extinction parameter, we had to take into account 
the (z - H) color variation, as a function of both age and 
of E(B-V). The MA05 synthetic stellar population mod- 
els with constant SFR, indicate that the (z - H) color is 
nearly constant in the redshift interval 1.4 < z < 2.5 at a 
given age, and E(B-V). In Figure [T5] we consider both 
the case of the same dust attenuation"ior H II regions and 
stellar continuum (EW(Ha), left panel), and th e case of 
extra attenuation, according the Equation [T] from Calzetti 
|eTaL] ( [2000]) (EW 00 (Ha), right panel, ct Section]?^ 
For comparison, we show the region occupied by syn- 
thetic stellar population models with constant SFR, so- 
lar metallicity, and five different values of dust extinc- 
tion (£(£- V)=0.0, 0.15, 0.2, 0.3, 0.4, chosen to include 
all values derived for the target galaxies). All the curves 
show the EW(Ha) and (z-H) color evolution in the age 
range 0.1-1.5 Gyr, and in steps of 0.1 Gyr (marked by 
small dots on the curves). The rest-frame EW(Ha) of a 
template at a given age and E(B - V) was computed as 
the ratio of the Ha flux, derived from the model SFR via 
the Kennicutt (1998) conversion, to the r-band flux den- 
sity of the templa te, assuming CSFR. For consistency, the 
Kennicutt (1998 ) conversion was corrected, as in Eq. (2), 
for the ChabrTer IMF, and the BC03 for the Chabrier IMF 
were used for the r-band continuum and the z-H colors. 

If one assumes that stellar continuum and nebular re- 
gions are equally extincted by dust, the equivalent width 
does not depend on dust attenuation. In fact, in the 
left panel of Figure [15] at a given age, and for increas- 
ing values of E(B-V), the model curves have redder 
colors, but the same EW(Ha) values. However, under 
this assumption, we do not find a good agreement be- 
tween data points and synthetic galaxy models. The mea- 
sured EW(Ha) tend to lie below the model curves with 
E(B-V) > 0, in regions that would correspond to un- 
physical ages older than that of the universe at z ~ 2, and 
null (i.e., black curve) or negative color excesses. On the 
other hand, if one assumes extra dust attenuation for neb- 
ular regions with respect to the stars (A Vi hii = Ay/0.44; 
right panel of Figure [15]) all the model curves tighten on 
a locus that forms an extension of the unreddened curve, 
and at a given age, £Wbo(Ha) exponentially decreases 
(as 10-°- 33 ^(i/o.44-i) _ 10 -i.72?(2MO) for increasing val- 
ues of E(B-V). From Figure 15 it is evident that, al- 
though with a large scatter, for zC-SINF galaxies the case 
of extra dust attenuation better matches with the model 
predictions. So, these results confirm what we found in 
the previous sections by comparing Ha luminosity with 
SFR from SED fitting. That is, for our sample Ha seems 
to be more extincted than the UV stellar continuum, in 



agreement with the prescription of |Calzetti et al. (2000). 
The fact that a substantial fraction of the data points in 
the right panel of Figure [15] show larger equivalent width 
than the model curves suggests that some of our galax- 
ies could have sub- solar metallicities, which would pro- 
duce higher equivalent widths (higher Ha luminosity) for 
a given SFR, age, and E(B-V) (cf.|Erb et al.|2006b|). 



5.1.3. Comparison of the integrated Ha properties of the 
zC-SINF sample with other IFS samples at z ~ 2 

Figures [T6}|T8] show the distributions of Ha properties 
of the zC-SINF sample (red filled circles). The integrated 
velocity dispersions cr(Ha)i nte grated and the half-light radii 
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Fig. 15. — The two panels show (z-H) colors vs Ha equivalent widths (EW in logarithmic scale) for the zC-SINF sample (red 
filled circles), by assuming that nebular regions suffer the same (left panel) or ^double dust extinction (right panel) with respect to 
the stars which generate the continuum. The reported uncertainties correspond to the formal measurement uncertainties, and do not 
take into account other sources of uncertainties (such as those from the absolute flux calibration of the SINFONI and broad-band 
photometry, or from aperture matching). Equivalent widths of Ha undetected sources are plotted at their 3 a limits (cf. Table|6}. The 
curves show predicted EW and colors for BC03 stellar population models with constant SFR and Chabrier IMF, and different values 
of dust extinction (E(B-V)), as labeled in the left panel. Dots on the curves represent increasing ages in steps of 0.1 Gyr, in the range 
0.1-1.5 Gyr. 
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Fig. 16. — Left panel: the integrated velocity dispersion a(Ka)i ntegr ated vs. the Ha flux F(Ha), where a (Ha) integrated is defined as the 
FWHM of the single Gaussian best fitting the integrated Ha profile of individual galaxies, and is corrected for instrumental spectral 
resolution. Besides galaxies in the zC-SINF sample (red filled circles) also galaxies in the IFS samples of FS09 (i.e., the SINS sample, 
open squares), and of L09 (green triangles) are displayed. The three histograms (arbitrarily normalized) show the relative distributions 
in F(Ha) and a(Ha) integra ted for the three samples, as labeled on the panel, and exclude the upper limits. Right panel: the half-light 
radius, r\j 2 , of the same galaxies vs. F(Ha), where r\/ 2 is measured as the radius of the circular aperture centered on the centroid of 
the Ha emission map, including 1/2 of the Ha luminosity, and determined based on the curve-of-growth. Symbols and histograms 
are the same as in the left panel. In the two panels the dashed lines mark out the line widths (left panel) and the sizes (right panel) 
above which the galaxies would be undetected (i.e., S/N < 3 per spectral or spatial resolution element, respectively) in our data set, 
as a function of Ha flux. 
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ri/2(Ha) are plotted as a function of the total observed 
fluxes F(Ha), of the luminosities L(Ha) (uncorrected for 
dust extinction), and of the stellar masses M*. Stellar 
masses are from BC03+CSFR models and Chabrier IMF, 
to allow a better comparison with other IFS samples us- 
ing these assumptions (see below). We also indicate in 
the Figures the typical sensitivity limits in terms of size 
and velocity width. These limits were determined follow- 
ing exactly the same procedure as described in FS09. In 
brief, we determined the average S/N per spatial resolu- 
tion element within the half-light radius for each galaxy 
based on the Ha maps. We then calculated the neces- 
sary increase in half-light radius for the S/N per resolu- 
tion element to drop below S/N = 3, keeping the total Ha 
flux and integrated line width fixed. These derived sizes 
correspond to the actual surface brightness sensitivity of 
each data set. The average limit in terms of r 1 / 2 (Ha) as 
a function of F(Ha) was obtained by taking the running 
median through the individual limits (and similarly as a 
function of L(Ha) and M*). The average limits in terms 
of cr(Ha)i ntegr ated were computed in a similar manner, cal- 
culating the line width at which the S/N per spectral res- 
olution element drops below S/N = 3 at fixed F(Ha) 
and r 1 / 2 (Ha). Because the limits show significant scatter 
among the objects, and some of the sources have an aver- 
age S/N per resolution element somewhat below 3, some 
of the individual measurements lie above the dashed lines 
in Figures 16 to 18 which represent typical sensitivity 
limits. 

The zC-SINF data are compared with the global Ha 
properties of other galaxy samples in the range 1.4 < 
z < 2.5 also observed with integral field spectrographs, 
taken from the SINS survey with SINFONI (FS09, open 
squares) and the samp le observed with OSIRIS at the 
Keck II telescope by |Law et al.| ( |2009| hereafter L09; 
green filled triangles). The targets for all three studies 
were drawn from parent samples at z ~ 2 selected with 
photometric criteria designed to identify massive star- 
forming galaxies, and are subject to similar biases intro- 
duced by the need for an accurate optical spectroscopic 
redshift and a sufficiently bright Ha emission. Although 
we include the results from L09, one should keep in mind 
that all targets of this sample were observed with the aid 
of AO and this may introduce differences in the derived 
Ha properties (see below). The median properties for 
these three samples (obtained by considering only the 
Ha detected objects) are shown in Table [7] The com- 
parisons show that the zC-SINF and SINS samples span 
nearly the same ranges in integrated Ha flux, velocity 
width, and half-light radius. They also overlap largely 
in terms of total Ha luminosity and stellar mass. Instead, 
although the L09 sample covers roughly the same range 
in Ha flux and luminosity, and in stellar mass, the in- 
tegrated velocity widths and half-light radii are system- 
atically and significantly smaller, by factors of approxi- 
mately two for the ensemble of objects, as well as for a 
given Ha flux/luminosity or a given stellar mass. As dis- 
cussed by FS09 and L09, these differences could arise 
from distinct population properties, in particular, more 
compact sizes. We emphasize that no selection based on 
sizes was applied for our zC-SINF sample and for the ma- 
jority of the SINS sample. The differences could also be 
partly due to a combination of the higher resolution and 
poorer surface brightness sensitivity for the AO-assisted 



OSIRIS data. Indeed, this is quite likely the case, given 
that smaller sizes for given stellar mass would have im- 
plied higher velocities, which instead are lower (see Fig- 
ure 18). This argues for the outer, high-velocity but low 
surface brightness regions having been lost. Moreover, 
the very similar M* distributions between the zC-SINF 
and L09 samples suggest that if the differences in sizes 
and velocity dispersions are attributed to intrinsic galaxy 
properties, stellar mass may not be the most important 
parameter driving these differences. 

The new zC-SINF data follow a similar trend of in- 
creasing integrated Ha velocity dispersion with increas- 
ing Ha flux, Ha luminosity, and stellar mass as seen 
among the SINS and OSIRIS sample galaxies. Like- 
wise, the Ha half-light radii tend to be larger for brighter 
and more massive sources, although the scatter is more 
significant. As discussed by FS09, the upper envelope 
of the distributions of cr(Ha)i nte grated and r^CHa) versus 
F(Ha), L(Ha), and likely results from observational 
biases. More specifically, for SINFONI seeing-limited 
observations, the detection limits for 1 hr integration lies 
just above the locus of the measurements, implying that 
the zC-SINF sample may be missing the most extended 
sources and/or those with broadest velocity widths at a 
given Ha flux/luminosity or stellar mass. However, such 
a bias cannot explain the lack of sources with small ve- 
locity dispersion at high Ha fluxes/luminosities and stel- 



lar masses seen in Figs. [16] T7J and[T8j which could be 
the consequence of a real/physical cr(Ha)i nte g r ated _ ^ re- 
lation. This issue, already pointed out by FS09 for the 
SINS sample, may not be surprising given that both kine- 
matic components that provide the total dynamical sup- 
port of galaxies (random motions and rotational/orbital 
motions) contribute to the measured cr(Ha)i nte grated, which 
thus probes the full gravitational potential of the galaxies. 
Besides being due to inclination effects, the scatter in the 
cr(Ha)i n tegrated _ diagram may possibly reflect varia- 
tions in the gas and dark matter mass fraction among the 
sources. 

6. KINEMATIC CLASSIFICATION 

For the zC-SINF pre-imaging data, the S/N ratio is gen- 
erally too low for a reliable quantitative kinematic classi- 
fication of the objects based on velocity fields and ve- 
locity dispersion maps. The distinction betwe en disk- 
like and m erger-like systems through kinemetry ( Shapiro 
et al.|2008] ) and detailed kinematic modeling will be per- 
formed based on the higher resolution data taken as part 
of the follow-up AO observations of the program. The 
results for the first few objects observed with AO in early 
stages of our on-going SINFONI LP and previous pro- 
grams are presented elsewhere (Shapiro et al. 2008; Gen- 
zel et al. 2008, 2011; Cresci et al. 2009). Here, we 
attempt a first crude kinematic classification in "rotation- 
dominated" and "dispersion-dominated" galaxies. 

Ideally, the distinction between rotation- and 
dispersion-dominated objects would rely on the ra- 
tio of the inclination-corrected circular or orbital velocity 
and the local intrinsic velocity disp ersion, v c [op (with 
the transition at v c / <7q ~ 1; e.g., |Genzel et al. 2008 
|Cresci et al.| |2009| ). Alternatively, the comparison 
between the contributions to the dynamical mass from 
random and from rotational/orbital motions can also 
be used (e.g., Epina t et al.||2009"] ). Both methods, 
however, require modeling of high quality kinematic 
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maps. Other indicators have been used in the literature, 
notably based on averaging over individual pixels or 
resolution elements across the galaxies t o estimate the 
intrin sic velocity dispersion (e.g., L09, |Green et al. 



2010| ). However, as demonstrated by |Davies et al.|20TT 
these estimates remain largely biased by the surface 
brightness sensitivity of the data and by beam-smearing 
effects. 

We therefore resorted to an alternative approach that 
can be followed for the majority of the zC-SINF galaxies 
with the current no- AO data sets. We used the working 
criterion introduced by FS09 and based on the ratio of 
half of the maximum observed velocity difference across 
the source (uncorrected for inclination) and the integrated 
velocity dispersion, v b S / 2 g(H af)i ntegr ated (each measured 
as described in Section |4.3.2 and reported in Table [5]). 
Based on simulations of disk galaxies with a range of in- 
clinations, sizes, and masses, as well as of beam smearing 
appropriate for the SINS sample, and thus also for the zC- 
SINF sample given the simil ar obse rving conditions and 
global Ha properties (Section |5.1.3| ), an intrinsic v c /<7o ~ 
1 corresponds roughly to v b S /2 0"(HaOintegrated ~ 0.4. This 
criterion was devised in the framework of rotating disks; 
since we do not distinguish between disk- and merger- 
like kinematics here, a v bs/2cr(Ha)i nt e g rated ^ 0.4 could 
also indicate dominant bulk orbital motions for merg- 
ing/interacting systems. 

Figure 
mass for 
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shows the v bs/2cr(Ha)i nt egrated versus stellar 
the subset of 24 zC-SINF objects for which we 
could derive v bs or an upper limit thereof. The measure- 
ments for 46 objects from the SINS Ha sample for which 
this ratio could be reliably determined are al so shown, 
along with those taken from the sample of |Law et al. 
2009). According to the criterion above, 13 of the 24 



objects have measured ratios or upper limits that lie in 
the rotation-dominated regime (including two within 1 a 
of the dividing line) and 11 objects fall in the dispersion- 
dominated regime (with three within 1 a of the dividing 
line). Bearing in mind that the v O b S /2cr(Ha0integrated ra- 
tio is a crude indicator, and that for many of the sources 
the velocity gradients can only be estimated over lim- 
ited regions with sufficiently high S/N, it is clear that the 
above results should be considered with the due caution. 
Nevertheless, the v bs/2cr(Ha)i ntegrate( j estimates suggest 
that roughly half of the zC-SINF galaxies appear to have 
kinematics dominated by disk rotation or bulk orbital mo- 
tions. We point out that these candidates also include two 
objects from the zCOSMOS "Pilot Sample," already clas- 
sified as rotating disks by means of kinemetry and kine- 
matic modeling (i.e., ZC407302, and ZC410542; Shapiro 
et al. 2008; Cresci et al. 2009) 3 . 

The relative proportions for the zC-SINF sample are 
roughly comparable to those inferred for the SINS Ha 
sample though with some small differences. Using the 
same v o bs/2cr(Ha0integrated indicator, about one-third of 
the SINS Ha sample has kinematics dominated by dis- 
persion, and two-thirds has kinematics dominated by 
rotational/orbital motions (including the kinematically- 
identified mergers). This is also reflected in the median 
v b S /2cr(Ha)i nt eg r ated of 0.44 (including upper limits) for 

3 For one of the other two Pilot Sample objects, ZC407376, the S/N 
of the seeing-limited SINFONI data is insufficient for quantitative anal- 
ysis through kinemetry and dynamical modeling, and the fourth one, 
ZC4101 16, is undetected in Ha (Table|5j. 
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Fig. 19. — The ratio of the maximum observed velocity dif- 
ference across the source, v bs, to twice the integrated velocity 
dispersion, a(Ha)i ntegra ted, as a function of galaxy stellar mass. 
The figure illustrates the working criterion of FS09, used here 
for a preliminary kinematic classification of the zC-SINF galax- 
ies. The Results for the subset of 24 sources for which the S/N 
allowed us to derive v bs (or an upper limit thereof) are com- 
pared with those from th e SINS Ha sample, and from the sam- 
ple of Law et al. (2009). Symbols are the same as in Figures 



16[ 17] and [T8] The dashed line (i.e., v bs/2cr(Ha)i n tegrated = 0.4) 
represents tneboundary between the dispersion-dominated and 
rotation-dominated regime. Galaxies above this line are tenta- 
tively identified as "rotation-dominated" according to the crite- 
rion developed by FS09. 

zC-SINF, and 0.57 for SINS (or 0.54 when excluding the 
five objects for which the ratio was measured from AO 
data). In contrast, the same indicator for the Law et al. 
(2007, 2009) sample observed with AO implies a me- 
dian of 0.28 and two-thirds of dispersion-dominated ob- 
jects. As discussed by FS09, these differences could be 
due to intrinsic differences between the samples studied 
and/or to sensitivity limitations; in particular, there may 
be a weak trend of increasing proportion of dispersion- 
dominated objects with decreasing stellar mass but the 
scatter of the data is too large to assess the trend robustly. 
While clearly more details about the kinematic proper- 
ties will be extracted from sensitive AO-assisted observa- 
tions and will allow an empirical calibration of this indi- 
cator, seeing-limited observations remain important and 
are less expensive in terms of observing time. In the fu- 
ture, the multi-IFU instrument KMOS at the VLT will 
provide much larger samples with kinematics but only 
in seeing-limited mode. The results from the zC-SINF 
no- AO data (and those from SINS) thus constitute valu- 
able references. The analysis above suggest a significant 
fraction of objects dominated by rotation/orbital motions, 
and it will be interesting to explore whether the kinematic 
mix varies as a function of global galaxy parameters once 
much larger samples allow to investigate this issue. 

7. CONCLUSIONS 

This paper presents the first results from the 'zCOS- 
MOS - SINFONI project', based on an ESO VLT Large 
Program aimed at studying the kinematical properties of 
z ~ 2 star-forming galaxies, as mapped by the Ha emis- 
sion line, with the near-infrared SINFONI AO-assisted 
integral field spectrograph. The zCOSMOS-SINFONI 
observing strategy includes (i) 'pre-imaging' natural see- 
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ing observations of 30 galaxies (i.e. the so-called zC- 
SINF sample), consisting in 1 hour of integration time 
per object, and (ii) deeper AO follow-up of ~ 20 sources, 
selected based on the pre-imaging results. 

The main distinct feature of this project with respect to 
the past consists in having assembled an unprecedented 
sample of galaxies at z ~ 2, for which high quality kine- 
matics resolved on spatial scales as small as ~ 1 kpc will 
be soon available. The zC-SINF sample selection in the 
zCOSMOS-Deep survey, offered the advantage of reli- 
able optical spectroscopic redshifts, and very wide area 
helping to cull targets with suitably bright and nearby 
stars enabling AO follow-up, as well as of the extensive 
multi-wavelength information from the COSMOS sur- 
vey. Moreover, our sample is uniform in terms of color 
criteria, quality of optical spectroscopic follow-up, and 
coverage in M* and SFR. Last, but not least, the 20 galax- 
ies observed with AO, will roughly double the current 
number of "benchmark" galaxies at z ~ 2 with spatially 
resolved kinematics on ~ 1 kpc scale. 

In this work the selection criteria used to assemble 
the full zC-SINF sample, the general galaxy proper- 
ties derived from SED fitting analysis, and the results 
from the first part of the SINFONI program, based on 
natural- seeing observations, were described. The results 
were then compared with the literature available for the 
broader population of star- forming galaxies at z ~ 2. The 
main conclusions are summarized in the following. 

1. The zC-SINF sample was culled from the zCOS- 
MOS survey, so to span a wide mass range (3 x 
10 9 < (M*/M©) < 2 x 10 11 ), to give the best possi- 
ble representation of the star-forming galaxy pop- 
ulation at z ~ 2. We showed that the sample is 
inevitably biased by the zCOSMOS pre-selection 
at Bab ^ 25, needed to obtain an adequate S/N in 
the VIMOS spectra, and by the lower limit that 
we imposed to the SFR (~ 1OM yr~ l ) of the se- 
lected objects, to match the typical detection limit 
for Ha line emission of z ~ 2 galaxies in 1 hour of 
pre-imaging observations. However, the sample in- 
cludes galaxies in a fairly wide range of SFRs (i.e., 
-40-300 M /yr). 



2. Using the |McCracken et al.|d20T0| ) i^-band selected 
multi-wavelength catalog, for each object we de- 
rived stellar mass, SFR, dust extinction Ay, and 
age, by fitting the photometric SED. We checked 
the results against different libraries of synthetic 
stellar population models, and different star for- 
mation histories (i.e., constant SFR, and "inverted- 
t" models), and found that zC-SINF masses, and 
SFRs, always agree within ~ 0.2-0.3 dex. How- 
ever, SFRs decrease and stellar masses increase 
when switching from constant SFR to inverted-r 
models, in such a way that the latter give a spe- 
cific SFR ~ 3 times lower than the former models. 
For various astrophysical reasons, and for giving 
SFRs in better agreement with those derived from 
the Ha luminosity, for the z ~ 2 galaxies in our 
sample we favor inverted-r models. In any way, 
secularly increasing SFRs are to be preferred com- 
pared to constant or declining SFRs. We used the 
Ilbert et al. ( 2009) /-band selected catalog to check 
also the SED fitting results against different photo- 



metric catalogs, and also found that they are stable 
within — 0.2 dex. 

3. Out of the initial 30, we detected 25 galaxies with 
SINFONI, and extracted Ha flux, integrated veloc- 
ity dispersion, and half-light radius (i.e., F(Ha), 
cr(Ha)i ntegrate d, and r^Ha)). The Ha integrated 
properties were then compared with those derived 
from IFS for other samples at z ~ 2. We found 
that the zC-SINF and the SINS sample (FS09; also 
based on SINFONI data) span the same range in 
F(Ha), cr(Ha)i n tegrated, and r!/ 2 (Ha). The fact that 
the results from this more uniform zC-SINF sam- 
ple are similar to those found with the more het- 
erogeneous SINS Ha sample provides an impor- 
tant confirmation in terms of the general properties 
of massive star- forming galaxies in the relevant M* 
and SFR range. 

4. Based on the Ha integrated fluxes, Ha luminosities 
(L(Ha)) were derived, and used to determine SFRs 



through the |Kennicutt| ( |1998| ) relation. SFRs were 
also derived from the UV luminosity (L1500), red- 
shifted in the 5-band for galaxies at z ~ 2, as well 
as from the SED fitting, and then compared with 
those from Ha. In this comparison, we considered 
two possibilities to account for dust extinction: (i) 
the HII regions are as extincted as the stellar con- 
tinuum; (ii) they a re ~ 2 times more ex tincted, as 
parametrized by Calzetti e t al.| ( |2000j ). For zC- 
SINF galaxies, we found that the additional dust 
extinction for nebular regions produces a better 
agreement among the Ha-derived SFR and other 
indicators. 

These results from the comparison of SFR(Ha) 
with other SFR estimates are further supported by 
our analysis of the Ha equivalent widths. Surface 
brightness limitations may affect our conclusions, 
but we estimate that these limitations are unlikely 
to account for the entire systematic offset towards 
lower SFR(Ha)s and EW(Ha)s that emerge when 
assuming A ViHII = A ViS ed- The evidence from our 
zC-SINF sample for extra attenuation towards the 
HII regions with respect to the bulk of the stel- 
lar population dominating the broad-band contin- 
uum emission agrees with earlier findings follow- 
ing similar arguments for the SINS Ha sample 
(FS09), and with the first direct measurements of 
the B aimer decrement in small samples of other 
massive, actively star-forming galaxies at z ~ 2. 

5. We presented a preliminary kinematic classifica- 
tion using a working criterion developed by FS09, 
and based on the ratio of two observed quantities, 
i.e., half of the maximum observed velocity dif- 
ference across the source (uncorrected for incli- 
nation) (v b s ), and the Ha line width (i.e., the in- 
tegrated velocity dispersion, cr(Ha)i nte grated)- Fol- 
lowing this criterion, we found that, out of the 24 
galaxies for which S/N ratio allowed us to retrieve 
a velocity gradient value, 14 objects appear to be 
rotation-dominated disk candidates, comparable to 
the fraction inferred in other IFS samples at high-z. 
The higher S/N and higher spatial resolution of the 
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follow-up AO-assisted SINFONI data, to be pre- 
sented in upcoming papers, will allow a more ro- 
bust classification and a detailed assessment of the 
reliability of our pre-classification. 
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APPENDIX 
APPENDIX A 

SPECTRAL ENERGY DISTRIBUTIONS FOR THE ZC-SINF SAMPLE 

Figure |20] shows the best-fit spectral energy distributions (SED) for the 30 galaxies of the zC-SINF sample. The 
open circles represent the photometric data points from to IRAC-bands. In each panel the best- fit curves obtained 
from BC03 constant SFR models are shown in black, together with those from MA05 constant SFR (green), and MA05 
inverted-r (red) models. The best- fit curves were obtained by fixing the redshift to the spectroscopic value, as measured 
from the Ha emission line (i.e., z(Ha)) for the 25 Ha-detected objects. For the remaining 5 Ha undetected objects the 
optical zCOSMOS redshift was used. Best fit parameters are shown in Table [4] 

APPENDIX B 

Ha MAPS, POSITION-VELOCITY DIAGRAMS AND INTEGRATED SPECTRA 

Figure [21] shows the /-band emission from HST+ACS (F814W), the velocity-integrated Ha linemap, the position- 
velocity diagram, and the integrated spectrum for each of the 25 Ha-detected galaxies of the zC-SINF sample. The 
Ha maps, position-velocity diagrams, and spectra were derived from 1 hour integration time for all objects, except 
ZC405081 (1.7 hour), and ZC407302 and ZC407376 (2 hours), and natural-seeing SINFONI observations (no-AO). 

The position-velocity (p-v) diagrams were extracted from the data cubes, without additional smoothing from median- 
filtering, using a synthetic slit 6 pixels wide (0"75) oriented along the major axis of the galaxies, and indicated by the 
rectangles on the Ha maps. Whenever possible, the galaxy major axis was identified based on the velocity field, as the 
direction of the steepest observed velocity gradient across the source (Ha kinematic major axis). For those sources with 
no clearly apparent velocity gradient, we took the morphological major axis. It is worth noticing that for the galaxies 
with a clear veloci ty gradient the kinematic major axis in general was found to be consistent with the morphological one 
(see Section |4.3.2| ). 

The integrated spectra were extracted from the unsmoothed data cubes in circular apertures shown (as black filled 
circles) on the Ha maps, using the aperture radii (r ap ) listed in Table [5] r ap were derived base d on the curve-of-growth 
method, so to enclose the total Ha flux (i.e., more than 90%), as also explained in Section [43 
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TABLE 1 

zCOSMOS-SINFONI 'pre-imaging' observing runs. 



Source 




RA 






DEC 




Date 


Period 


Band 


PSF FWHM [arcs] 


A int L a J 


ZC407376 a 


10: 


00 


:45 


.1 


02 


:07: 


:05 


.2 


16-23 Apr 2007 


P79 


K 


0.60 


7200 


ZC407302 a 


09: 


59 


:56 


.0 


02 


:06: 


:51 


.3 


16-23 Apr 2007 


P79 


K 


0.56 


7200 


ZC410116 a 


09: 


59 


:36 


.5 


02 


:16: 


:14. 


.5 


18-20 Aug 2007 


P79 


K 


0.71 


3600 


ZC410542 a 


09: 


59 


:54 


.7 


02 


:17: 


:48. 


.3 


18-20 Aug 2007 


P79 


K 


0.82 


3600 


ZC409985 


09: 


59 


:14 


.2 


02 


:15: 


:46. 


.9 


21-22 Feb 2008 


P81 


K 


0.69 


3600 


ZC410123 


10: 


02 


:06 


.5 


02 


:16: 


:15 


.6 


22-23 Feb 2008 


P81 


K 


0.29 


3600 


ZC405226 


10: 


02 


:19 


.5 


02 


:00: 


:18. 


.2 


22-23 Feb 2008 


P81 


K 


0.62 


3600 


ZC404073 


10: 


01 


:42 


.8 


01 


:56: 


:12 


.7 


22-23 Feb 2008 


P81 


K 


0.57 


3600 


ZC403027 


10: 


02 


:21 


.3 


01 


:52: 


:34. 


.3 


23-24 Feb 2008 


P81 


K 


0.70 


3600 


ZC400528 


09: 


59 


:47 


.6 


01 


:44: 


:19. 


.0 


26-27 Feb 2009 


P83 


K 


0.49 


3600 


ZC400569 


10: 


01 


:08 


.7 


01 


:44: 


:28. 


.3 


26-27 Feb 2009 


P83 


K 


0.60 


3600 


ZC403103 b 


09: 


59 


:19 


.0 


01 


:52: 


:50 


.1 


1-2 Mar 2009 


P83 


K 


0.64 


3600 


ZC405545 


10: 


00 


:32 


.8 


02 


:01: 


:17. 


.7 


14-15 Mar 2009 


P83 


K 


0.84 


3600 


ZC415876 


10: 


00 


:09 


.4 


02 


■36 


:58. 


.3 


14-15 Mar 2009 


P83 


K 


0.49 


3600 


ZC405081 


10: 


00 


:22 


.7 


01 


:59: 


:46. 


.2 


16-17 Mar 2009 


P83 


K 


0.54 


6000 


ZC401925 


10: 


01 


:01 


.7 


01 


:48: 


:38. 


.2 


16-17 Mar 2009 


P83 


K 


0.52 


3600 


ZC404221 


10: 


01 


:41 


.3 


01 


:56: 


:42. 


.8 


16-17 Mar 2009 


P83 


K 


0.55 


3600 


ZC403741 


10: 


00 


:18 


.4 


01 


:55: 


:08. 


.1 


20-21 Mar 2009 


P83 


H 


0.60 


3600 


ZC404987 


10: 


00 


:40 


.4 


01 


:59: 


:27. 


.7 


24-25 Mar 2009 


P83 


K 


0.58 


3600 


ZC412369 


10: 


01 


:46 


.9 


02 


:23: 


:24. 


.6 


24-25 Mar 2009 


P83 


K 


0.56 


3600 


ZC406690 


09: 


:58 


:59 


.1 


02 


:05: 


:04. 


.2 


29-30 Dec 2009 


P84 


K 


0.68 


3600 


ZC411737 


10: 


00 


:32 


.4 


02 


:21 


20. 


.9 


01-02 Jan 2010 


P84 


K 


0.48 


3600 


ZC4 13597 


09: 


59 


:36 


.4 


02 


:27: 


:59 


.1 


01-02 Jan 2010 


P84 


K 


0.54 


3600 


ZC405501 


09: 


59 


:53 


.7 


02 


:01: 


:08. 


.8 


25-26 Jan 2010 


P84 


K 


0.43 


3600 


ZC407928 


09: 


59 


:53 


.4 


02 


:08: 


:41. 


.5 


25-26 Jan 2010 


P84 


K 


0.47 


3600 


ZC415087 


10: 


00 


:32 


.4 


02 


:33: 


AO. 


.6 


08-09 Feb 2010 


P84 


K 


0.76 


3600 


ZC4 13507 


10: 


00 


:24 


.2 


02 


:27: 


■41. 


.3 


10-11 Feb 2010 


P84 


K 


0.43 


3600 


ZC409129 


10: 


01 


:28 


.0 


02 


:12: 


A6. 


.6 


07-08 Mar 2010 


P84 


K 


0.80 


3600 


ZC410041 


10: 


00 


:44 


.3 


02 


:15: 


:58. 


.4 


07-08 Mar 2010 


P84 


K 


0.78 


3600 


ZC405430 


10: 


01 


:32 


.9 


02 


:00: 


:56 


.5 


01-02 Apr 2010 


P84 


K 


0.65 


3600 



a Targets observed as a 'pilot sample' of the zCOSMOS and SINS teams collaboration, that are also included in 
the SINS sample presented in FS09. They are named here with the new zCOSMOS IDs. They had provisional 
IDs in 2007, respectively, ZC-772759, ZC-782941, ZC-946803, and ZC-1 101592, and so they were named in 
FS09. 

b This galaxy is the only one in the zC-SINF sample being preselected with the UGR/BX criterion in 
zCOSMOS-Deep, instead of the BzK criterion (see also L07, Sect.[33] and Fig. [4}. 
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TABLE 2 



Optical photometry of the zC-SINF sample 


ID 


u* 


Bj g + Vj r + 


i + 





ZC400528 24.14 ±0.02 23.71 ± 0.01 23.64 ± 0.02 23.34 ± 0.01 23.13 ± 0.01 22.86 ± 0.01 22.80 ± 0.03 

ZC400569 24.84 ±0.05 24.27 ± 0.02 24.17 ± 0.03 23.79 ± 0.02 23.60 ± 0.01 23.28 ± 0.02 23.05 ± 0.04 

ZC401925 24.32 ±0.02 24.10 ±0.01 24.08 ± 0.02 23.72 ± 0.02 23.67 ± 0.01 23.63 ± 0.02 23.57 ± 0.05 

ZC403027 25.58 ±0.07 24.93 ± 0.03 24.82 ± 0.06 24.15 ± 0.03 24.09 ± 0.02 23.83 ± 0.03 23.55 ± 0.05 

ZC403103 24.97 ±0.04 24.28 ± 0.01 24.23 ± 0.03 23.72 ± 0.02 23.64 ± 0.01 23.40 ± 0.02 23.47 ± 0.05 

ZC403741 24.27 ±0.02 24.06 ± 0.01 24.12 ± 0.02 23.77 ± 0.02 23.53 ± 0.01 23.08 ± 0.01 22.65 ± 0.02 

ZC404073 25.95 ±0.10 24.99 ± 0.03 24.83 ± 0.05 24.12 ± 0.02 23.95 ± 0.02 23.71 ± 0.02 23.54 ± 0.06 

ZC404221 24.09 ± 0.02 23.86 ± 0.01 23.75 ± 0.02 23.49 ± 0.02 23.49 ± 0.01 23.35 ± 0.02 23.51 ± 0.06 

ZC404987 24.85 ± 0.03 24.57 ± 0.02 24.37 ± 0.02 24.06 ± 0.02 24.08 ± 0.02 23.97 ± 0.03 23.86 ± 0.07 

ZC405081 24.78 ± 0.02 24.44 ± 0.02 24.46 ± 0.03 24.06 ± 0.02 24.06 ± 0.02 23.91 ± 0.03 23.96 ± 0.08 

ZC405226 24.77 ± 0.03 24.21 ± 0.01 24.04 ± 0.03 23.58 ± 0.02 23.55 ± 0.01 23.33 ± 0.02 23.29 ± 0.04 

ZC405430 25.99 ±0.10 25.33 ± 0.04 25.18 ± 0.07 24.76 ± 0.04 24.69 ± 0.04 24.51 ± 0.05 24.37 ±0.12 

ZC405501 24.22 ± 0.02 23.97 ± 0.01 24.00 ± 0.03 23.61 ± 0.02 23.62 ± 0.02 23.41 ± 0.02 23.40 ± 0.06 

ZC405545 25.65 ±0.06 25.02 ± 0.03 24.95 ± 0.05 24.47 ± 0.03 24.50 ± 0.03 24.19 ±0.04 24.21 ±0.11 

ZC406690 22.98 ±0.01 22.58 ± 0.01 22.62 ± 0.01 22.15 ± 0.01 22.05 ± 0.00 21.88 ±0.01 21.89 ±0.02 

ZC407302 24.03 ± 0.01 23.62 ± 0.01 23.51 ± 0.01 23.05 ± 0.01 23.06 ± 0.01 22.85 ± 0.01 22.71 ± 0.02 

ZC407376 25.28 ±0.04 24.75 ± 0.03 24.73 ± 0.05 24.10 ±0.03 24.11 ± 0.03 23.78 ± 0.03 23.59 ± 0.07 

ZC407928 25.47 ±0.05 24.81 ± 0.02 24.75 ± 0.06 24.17 ± 0.03 24.09 ± 0.03 23.80 ± 0.03 23.82 ± 0.08 

ZC409129 25.66 ±0.06 24.85 ± 0.02 24.77 ± 0.04 24.32 ± 0.02 24.17 ± 0.02 23.87 ± 0.03 23.67 ± 0.06 

ZC409985 24.41 ±0.02 24.16 ± 0.02 24.29 ± 0.03 23.66 ± 0.02 23.55 ± 0.01 23.46 ± 0.02 23.49 ± 0.05 

ZC410041 24.71 ± 0.02 24.24 ± 0.02 24.24 ± 0.03 23.83 ± 0.02 23.82 ± 0.02 23.67 ± 0.02 23.70 ± 0.07 

ZC410116 25.41 ±0.05 25.31 ± 0.04 25.24 ± 0.06 24.69 ± 0.03 24.59 ± 0.03 24.15 ±0.03 23.81 ± 0.06 

ZC410123 24.62 ± 0.03 24.21 ± 0.01 24.23 ± 0.03 23.75 ± 0.02 23.84 ± 0.02 23.63 ± 0.02 23.61 ± 0.06 

ZC410542 23.97 ± 0.01 23.78 ± 0.01 23.82 ± 0.02 23.46 ± 0.02 23.26 ± 0.01 22.78 ± 0.01 22.32 ± 0.02 

ZC41 1737 24.70 ± 0.03 24.25 ± 0.02 24.25 ± 0.03 23.84 ± 0.02 23.77 ± 0.02 23.63 ± 0.02 23.74 ± 0.07 

ZC412369 24.37 ± 0.03 24.01 ± 0.01 23.89 ± 0.03 23.38 ± 0.02 23.40 ± 0.01 23.06 ± 0.02 22.88 ± 0.04 

ZC413507 25.64 ± 0.08 24.73 ± 0.03 24.78 ± 0.05 24.06 ± 0.03 24.02 ± 0.02 23.82 ± 0.03 23.71 ± 0.07 

ZC413597 25.28 ±0.06 24.82 ± 0.02 24.70 ± 0.04 24.18 ± 0.02 23.95 ± 0.02 23.90 ± 0.03 23.80 ± 0.07 

ZC415087 24.80 ± 0.03 24.56 ± 0.02 24.50 ± 0.04 24.04 ± 0.03 24.21 ± 0.03 23.86 ± 0.03 23.78 ± 0.08 

ZC415876 25.23 ± 0.04 24.64 ± 0.02 24.57 ± 0.04 23.99 ± 0.02 23.95 ± 0.02 23.72 ± 0.02 23.69 ± 0.06 
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TABLE 3 

NEAR-IR/MID-IR photometry of the zC-SINF sample 



ID 






J 








H 






K 




m 3 .6 




^4.5 




^5.8 


^8.0 


S24 [^Jy] 


ZC400528 


22 


.35 


± 


0. 


,11 


21 


.75 ± 


0.04 


21, 


,08 


±0.03 


20 


.43zb 0. 


10 


20.30 zbO.10 


20.21 zb 0.10 


20.29 zbO.15 


228 zb 18 


ZC400569 


22 


.11 


± 


0. 


.09 


21 


.29 zb 


0.03 


20. 


.69 


zb 0.02 


20 


.25zb 0. 


10 


20.05 zbO.10 


20.11 zb 0.11 


20.62 zbO.16 


141 zb 29 


ZC401925 


23 


.01 


± 


0. 


,14 


22. 


.92 ± 


0.12 


22, 


,74 


zb 0.12 


22. 


.52zb 0. 


11 


22.64 ±0.12 


22.11 zb 0.23 


22.99 zbl.06 


< 80 c 


ZC403027 


23 


.24 


± 


0. 


,15 


22. 


.39 ± 


0.05 


22, 


1 Q 


zz yj.yj i 


21 


.53zb0. 


10 


21.56 zbO.10 


21.53 zb 0.16 


21.76 zbO.35 


84 zb 16 


ZC403103 


22 


.91 


± 


0. 


13 


23. 


.94 ± 


0.20 


23, 


/IO 

Ay 


i n oq 
zt U.Z3 


23. 


.42zb 0. 


14 


23.87 zbO.25 


>21.3 b 


>21.0 b 


< 80 c 


ZC403741 


21 


.81 


± 


0. 


,05 


21 


.30 ± 


0.02 


21. 


no 
,UZ 


i r\ no 
zt U.UZ 


20 


.51zb0. 


10 


20.49 zbO.ll 


20.54 zb 0.11 


21.02zb 0.15 


< 80 c 


ZC404073 


22 


.63 


± 


0. 


,12 


22. 


.35 ± 


0.04 


21, 


,yo 


zt U.UD 


21 


.34zb 0. 


10 


21.20 zbO.10 


20.87 zb 0.12 


21.39zb 0.27 


82 zb 18 


ZC404221 


22 


.93 


± 


0. 


,17 


22. 


.83 ± 


0.08 


22. 


44 


ziz yj.yjy 


21 


.86zb 0. 


10 


21.80 zb 0.11 


21.78 zb 0.17 


21.03zb 0.19 


< 80 c 


ZC404987 


23 


.31 


± 


0. 


,15 


23. 


.19 ± 


0.11 


22. 


,77 


zb 0.13 


22. 


.27zb 0. 


11 


22.39 zb 0.12 


23.04 zb 0.52 


22.81zb 1.01 


< 80 c 


ZC405081 


23 


.27 


± 


0. 


,15 


23. 


.40 zb 


0.15 


22. 


.91 


zb 0.16 




-1.00 a 




-1.00 a 




-1.00 a 


-1.00 a 


< 80 c 


ZC405226 


23 


.57 


zb 


0. 


,44 


22. 


.61 ± 


0.07 


22, 


.33 


zb 0.09 


21 


.97zb 0. 


10 


22.07 zb 0.11 


22.32 zb 0.27 


21.73 zb 0.36 


< 80 c 


ZC405430 


25 


.04 


± 


1. 


,01 


23. 


.61 ± 


0.16 


23, 


.10 


zb 0.17 


23. 


.41zb0. 


14 


23.79 zb 0.25 


>21.3 b 


>21.0 b 


< 80 c 


ZC405501 


22 


.81 


zb 


0. 


,19 


23. 


.12 ± 


0.14 


22, 


9S 

,Z J 


-1-0 11 


22. 


.05zb 0. 


10 


21.99 zbO. 


11 


22.70 zb 0.39 


>21.0 b 


< 80 c 


ZC405545 


23 


.74 


± 


0. 


,26 


22. 


.99 zb 


0.12 


22, 


QO 


zt U.l 1 


22. 


.64zb 0. 


11 


22.50 zb 0. 


12 


23.65 zb 0.89 


>21.0 b 


< 80 c 


ZC406690 


21 


.41 


± 


0. 


,03 


21 


.11 ± 


0.02 


20, 


,o 1 


zz yj.yjD 


20 


.93zb 0. 


10 


20.89 zb 0. 


10 


20.72 zb 0.12 


20.84 zb 0.18 


209 zb 16 


ZC407302 


22 


.17 


± 


0. 


,05 


21 


.98 zb 


0.03 


21. 


,48 


zb 0.04 




-1.00 a 




-1.00 a 




-1.00 a 


-1.00 a 


295 zb 16 


ZC407376 


22 


.67 


± 


0. 


,10 


22. 


.32 ± 


0.05 


21, 


,79 


zb 0.06 


21 


.86zb 0. 


10 


21.69 zbO. 


10 


20.97 zb 0.12 


21.53 zb 0.27 


103 zb 16 


ZC407928 


23 


.62 


zb 


0. 


,23 


23. 


.25 ± 


0.12 


22. 


,63 
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a Sources with IRAC fluxes equal to -1.00 are blended in all the IRAC bands. This was deduced based both on the visual 
inspection of IRAC image stamps, and on the empirical criterion, tested on extensive simulations by the GOODS team, which 
consists in classifying as blended a ll the galaxies for which the angular separation between the IRAC and the ^T-band position 
exceeds 0"6 (cf. Da ddi et al.|2007) >. 

b For IRAC and undetected sources we list the 5a magnitude lower limits (Sand ers et al.|2 007 ), corresponding to the flux 
upper limits shown in Fig. |20| ( black arrows) 

c 5a detection limit for the JvTIPS 24 fim catalog ( |Le Floc'h et al.|2 009 ). 
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TABLE 7 

Median properties for the zC-SINF, SINS, and L09 samples compared in Figures [T6HT81 



Sample 


^*,med 

[1O 1O M ] 


F(Ho;) med 
[l(T 17 erg cm- 2 s" 1 ] 


L(Hai) m ed 

[lO^ergs- 1 ] 


<7 (H Q! ) integrated , med 

[km s _1 ] 


r 1/2 (HaO m ed 
[kpc] 


zC-SINF 


1.3 


7.9 


2.4 


123 


2.7 


SINS 


2.6 


11.0 


3.5 


130 


3.1 


L09 


1.2 


9.0 


3.4 


70 


1.4 
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Fig. 20. — Best-fit SEDs for the whole zC-SINF sample in units of Fa [erg/cm 2 /s/ A] vs observed wavelengths (A obs)[/xm]. Open cir- 
cles are the observed fluxes from the photometric data. The black, green, and red curves are the best-fit SEDs built with BC03+CSFR, 
MA05+CSFR, and MA05+inv-r SFH models, respectively. The black arrows are 2a flux upper limits (see also Table [3}. 
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Fig. 20.— Continue. 
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Fig. 21.— From left to right we show the /-band stamps from the COSMOS HST+ACS/F8 1 4W images, and the Ha line maps, 
position- velocity diagrams, and integrated spectra obtained from SINFONI natural seeing (no- AO) observations, for the Ha-detected 
zC-SINF sources. The circular aperture for the integrated spectrum and the synthetic slit used to extract the p-v diagrams along the 
major axis are indicated on the Ha maps. Vertical green hatched bars in the spectra show the locations of bright night sky lines, with 
width corresponding to the FWHM of the effective spectral resolution of the data. 
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